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SECTION  I 
INTRODUCTION 


The  work  done  under  Contract  No.  N00014-77-C-0296  is 
described  in  a  number  of  different  places  which  are  drawn  together 
and  summarized  in  this  final  report.  First,  the  publications  credited 
to  this  contract  which  are  available  from  the  open  literature  are 
listed  below. 

1.  "A  Score  for  Correct  Data  Association  in  Multi-Target  Track¬ 
ing."  by  D.  L.  Alspach  and  R.  N.  Lobbia.  appears  in  the 
December  1979  Proceedings  of  the  Decision  and  Control  Con¬ 
ference  held  in  Fort  Lauderdale.  Florida. 

2.  "Sound  Speed  Estimation  as  a  Means  of  Improving  Target 
Tracking  Performance."  by  D.  L.  Alspach,  G.  L.  Mohnkem, 
and  R.  N.  Lobbia.  appears  in  the  Proceedings  of  the  13th 
Asilomar  Conference  on  Circuits,  Systems,  and  Computers, 
November  1979,  Pacific  Grove,  California. 

3.  "Multiple  Coherence,"  by  R.  D.  Trueblood  and  D.  L.  Alspach, 
was  presented  and  appears  in  the  Proceedings  of  the  November 
1978  Conference  on  Systems  Science. 

4.  "A  Cfcse  Study  in  Adaptive  Sound  Speed  Estimation,"  by  R.  N. 
Lobbia  and  D.  L.  Alspach,  was  presented  at  and  appears  in 
the  Proceedings  of  the  November  1978  Asilomar  Conference  on 
Systems  Science. 

5.  "Data  Association  Algorithms  for  Large  Area  Surveillance,"  by 
C.  M.  Petersen  and  C.  L.  Morefield,  ORSA/TIMS  Meeting, 

May  1978. 

6.  "Application  of  0-1  Integer  Programs  to  Multi-Target  Tracking 
Problems,"  by  C.  L.  Morefield,  IEEE  Transactions  on  Automatic 
Control,  June  1977. 


Copies  of  the  papers  published  during  the  second  year  of  the 
contract  are  included  as  Volume  II  of  this  report.  Additional  copies 
of  all  the  above  have  been  transmitted  to  ONR. 

The  remainder  of  the  work  performed  under  this  contract  to 
ONR  is  of  two  different  types.  A  specific  approach  to  reducing  the 
computational  problem  in  implementing  a  certain  class  of  multi-target 
problems  is  contained  in  Section  II.  This  section  has  been  submitted 
in  essentially  this  form  to  Information  Sciences  for  publication. 

The  second  part  of  the  work  is  a  general  discussion  of  all 
aspects  of  the  overall  ocean  monitoring  and  control  problem.  This 
essentially  independent  white  paper  is  entitled  "A  Discussion  of  the 
Ocean  Tactical  Targeting  Problem"  and  is  included  as  Section  III  of 
this  final  report. 
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SECTION  II 

AN  EFFICIENT  TRACK  ASSOCIATION  ALGORITHM  FOR 
THE  MULTITARGET  TRACKING  PROBLEM 


ABSTRACT 

As  the  final  contribution  for  this  contract,  a  new  and  extremely 
promising  approach  is  presented  for  the  solution  of  the  multitarget 
tracking  problem  for  an  important  class  of  targets.  For  the  class  of 
targets  that  can  be  described  as  a  linear,  time-invariant  system,  it 
appears  that  this  approach  may  be  considerably  more  efficient,  compu¬ 
tationally.  than  existing  procedures. 

A  survey  of  approaches  to  the  solution  of  the  multitarget 
tracking  problem  has  been  presented  in  Reference  2.  As  indicated 
there,  most  algorithms  require  the  use  of  the  Kalman  filter  to  provide 
state  estimation  for  data  association  and  track  assignment.  In  this 
report,  the  need  for  state  estimation  to  identify  feasible  tracks  is 
eliminated.  The  use  of  simpler  input-output  models  for  track  pre¬ 
diction  results  in  a  very  substantial  reduction  in  computational 
burden  when  compared  with  methods  requiring  state  estimation. 

The  utility  of  the  method  presented  here  depends,  as  do 
Kalman  filter-based  methods,  upon  the  definition  of  data  windows 
that  permit  the  identification  of  feasible  tracks  and  the  elimination 
of  infeasible  tracks  in  a  timely  manner  and  according  to  statistically 
meaningful  criteria.  After  identifying  feasible  tracks,  the  Kalman 
filter  can  be  utilized  to  obtain  the  required  state  estimates  and  to 
analyze  and  refine  the  information  in  the  feasible  track  files. 
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1.0  INTRODUCTION  TO  THE  PROBLEM 


The  problem  of  tracking  multiple  targets  using  measurements 
whose  source  is  questionable  seems  to  be  Intrinsic  to  surveillance 
systems.  Problems  having  this  character  have  received  considerable 
attention  since  Sittler's  paper  [1]  in  1964.  A  recent  survey  paper 
[2]  by  Bar-shalom  has  described  the  problem  and  has  reviewed  the 
literature  and  principal  results  regarding  the  solution  of  the  pro¬ 
blem.  Our  discussion  will  assume  familiarity  with  this  body  of 
work  to  at  least  the  depth  provided  by  Reference  2. 

The  problem  that  is  considered  here  has  the  following  de¬ 
scription.  Measurements  are  obtained  at  discrete  times;  t^,  t^> 
•••»  t^,  .  At  each  measurement  time  a  collection  of  meas¬ 

urements 

-*k"  Ui(tk)j  tml*  2»  — *  “k  I*  k“i»  2»  ••• 

is  obtained.  Each  measurement  jZ^Ct^)  is  generated  by  a  different 
source.  For  this  discussion,  the  set  of  measurements  up  to  and  in¬ 
cluding  is  denoted  as 

i  -  |  *ii  * "  1»  2.  k  I 

No  ordering  of  measurements  is  preserved,  necessarily,  from  one 
sampling  time  to  another.  For  example,  _z^(t^)  and  js^tj) 
emanate  from  entirely  different  sources  with  the  result  that  they 
are  completely  independent  of  each  other.  It  is  this  lack  of 
ordering  that  produces  that  multitarget  tracking  problem  that  is 
to  be  addressed  here. 

If  there  were  no  uncertainty  regarding  the  source  of  each 
measurement,  the  multitarget  tracking  problem  would  be  solved  by 
a  straightforward  application  of  the  methods  developed  for  single¬ 
target  tracking  (e.g. ,  extended  Kalman  filter).  As  described  in 


» 


» 


I 


I 


I 


Reference  2,  previous  approaches  to  the  problem  have,  generally, 
been  based  on  state-space  models  and  recursive  filtering  methods. 

The  new  aspects  of  the  multitarget  tracking  problem  stem  from  the 
necessity  for  measurement-to-track  assignment  of  data  and  for  the 
incorporation  into  the  filtering  algorithm  of  the  possibility  of 
Incorrect  data  assignment.  It  is  to  this  track  assignment  problem 
that  attention  is  addressed  in  this  report. 

To  provide  a  basis  for  the  discussion  consider  the  following 
mathematical  model  for  the  problem.  Suppose  that  M^,  jk-1,  2,  •••} 
measurement  vectors  j^(t^) ,  {  i  •  1,  2,  •••,  M^}  are  obtained  at  meas¬ 
urement  time  t,  .  These  measurements  are  assumed  to  emanate  from 
k 

different  sources.  Thus,  there  are  sources  at  t^.  However,  only 
a  subset  of  these  sources  have  interest  to  the  survelllsnce  system 
although  the  interest  in  a  source  is  to  be  deduced  from  the  analysis 
of  the  data.  Briefly,  each  measurement  (t^)  may  be  either  a  sig¬ 
nal  ^(t^)  plus  noise  v^(t^)  or  noise  v^(t^) ,  alone.  The  solution 
of  this  detection  problem  is  fundamental  to  the  satisfactory  perform¬ 
ance  of  the  surveillance  system. 

The  detection  is  accomplished  through  the  analysis  of  data 
z.  obtained  over  some  interval  [t^,  t^].  Generally,  sources  of  in¬ 
terest  (l.e.ft  targets)  exhibit  motion  during  the  observation  interval. 
The  motion  is  characterized  by  defining  a  state-space  model  that  is 
based  on  a  differential  or  a  difference  equation.  For  this  analysis, 
the  target  is  assumed  to  be  described  by  a  state  vector  x,  satis¬ 
fying  a  linear  difference  equation 

iJ(tk+l)  "*2j(tk>  +I'^j(tk)  +£j(tk)  *  (1,1) 


The  vector  represents  a  control  input  to  the  system  which  is 
unknown  to  the  surveillance  system  but  cannot  be  neglected  (e*g*f 
heading  changes) •  The  vector  w^  represents  random  influences  on 
the  system  motion  and  is  assumed  to  be  described  as  a  white  noise 
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sequence  with  zero  mean  and  known  covariance  matrix  Q.  The  matrices 

Q,  r,  and  Q  are  assumed  to  be  known  and  time-invariant.  The  time- 
invariance  of  the  system  is  a  restriction  on  the  generality  of  the 
problem  described  by  Bar-shalom  [2]  in  his  survey  paper.  Nonetheless, 
it  does  describe  an  Important  class  of  problems  that,  based  on  [2], 
has  not  been  treated  extensively. 

The  i^  measurement  at  t^  is  related  to  the  target  state  in 
the  following  manner 

z^)  -  H  Xj (tk)  +  sfy}  <1.2) 

The  matrix  H  is  known  and  time- invariant.  The  noise  v.  is  assumed 

1 

to  be  zero  mean,  white  noise  with  time-invariant  covariance  matrix 

R.  Measurements  ^(t^)  that  are  not  generated  by  a  target  shall  be 
assumed  here  to  be  Independent  and  identically  distributed  random 
variables  drawn  from  a  uniform  distribution  defined  on  the  measure¬ 
ment  space. 

In  Equation  (1.2),  the  1^  measurement  z  (t.)  is  indicated  as 

£tl  *  * 

being  related  to  the  state  of  the  j  target.  If  the  specific  tar¬ 
get  were  known,  the  multitarget  tracking  problem  would  be  readily 
solved.  The  problem  to  be  addressed  here  relates  to  the  data  asso¬ 
ciation  problem  of  establishing  the  specific  target  j  to  which  the 
1th  measurement  relates. 

Several  approaches  have  been  taken  to  the  solution  of  the 
data  association  problem.  Basic  to  most  of  these  is  the  application 
of  Kalman  filter  for  state  estimation.  The  Kalman  filter  is  used 
to  estimate  the  state  at  t^  from  the  data .assigned  to  a  candidate 
track.  Then,  this  estimate  is  used  to  predict  the  next  measurement 
z^tj^).  Typically,  some  type  of  data  window  is  defined  to  establish 
the  next  measurement  that  is  to  be  assigned  to  the  track.  When  more 
than  one  measurement  could  be  assigned,  various. alternatives  arise 
ranging  from  accepting  only  a  single  "or  best"  measurement  to 


"splitting"  the  track  into  several  candidates.  Since  the  possibil¬ 
ity  exists  that  a  measurement  might  be  mistakenly  assigned  to  a 
track,  some  researchers  have  suggested  modifications  to  the  Kalman 
filter  to  reflect  the  possible  misassignment .  It  is  important  to 
recognize  that  several  different  Kalman  filters  may  operate  at  each 
time;  several  more  than  the  number  of  actual  target  tracks.  Thus, 
computational  burden  can  become  an  important  concern.  Several  re¬ 
searchers  have  explored  methods  for  establishing  "feasible"  tracks 
or  "most  likely"  tracks  as  a  screening  mechanism  to  eliminate  tracks 
to  reduce  the  computational  burden. 

In  the  next  section,  an  approach  to  the  problem  of  determin¬ 
ing  feasible  tracks  that  seem  to  be  extremely  efficient  for  the 
class  of  linear,  time-invariant  systems  defined  above.  In  brief, 
the  approach  is  based  on  the  recognition  that  input-output  models 
can  be  defined  that  eliminates  the  need  for  state  estimation  to 
accomplish  the  data  association  problem.  If  state  estimation  can 
be  deferred  until  after  establishing  feasible  tracks,  considerable 
unnecessary  filter  computation  can  be  eliminated. 


2.0  DATA  ASSOCIATION  USING  INPUT-OUTPUT  MODELS 

Consider  the  dynamic  and  measurement  model  for  a  specific 
target.  For  convenience,  the  subscripts  i  and  j  appearing  in  Equa¬ 
tions  (1.1)  and  (1.2)  are  omitted. 

x(k+l)  “  $x(k)  +  r u(k)  +  w(k)  (2.1) 

z(k)  -  Hx(k)  +  v(k)  (2.2) 

It  is  straightforward  to  eliminate  the  state  variables  and 
to  determine  a  difference  equation  model  involving  only  the  Inputs 
and  outputs  when  the  target  is  operating  in  steady-state.  By 
eliminating  the  state  variables,  the  need  for  estimating  them  is 
eliminated.  The  input-output  model  can  be  used  directly  to  pre¬ 
dict  measurements  at  future  times.  The  variance  of  the  error  in 
the  predictions  can  be  determined  and  used  to  establish  feasible 
tracks . 

To  establish  an  input-output  model,  let  q  denote  a  forward 
difference  operator. 

x(k+l)  ^  qx(k)  (2.3) 

Then,  (2.1)  can  be  written  as 

(ql-$)x(k)  -  Tu(k)  +  w(k) 

Solving  for  x(k) ,  one  obtains 

x(k)  -  (qI-$)-1tru(k)  +w(k)] 

Using  this  result  in  (2.2),  it  follows*  that 

*The  definition  A  *  ■ 


is  used. 


..J|UIJLWUMUW>II«UWI!^ 


z(k)  -  H  (ql-  ^  )  _1  [  r  u(k)  +  w(k)] 
or 

det(ql- $)£(k)  »  Hadj(ql~$)  [  Tii(k)  +  w(k)J 

+  det(ql- $)v(k)  (2.4) 

where  adj  (ql-  $)  denotes  the  adjoint  of  the  matrix  (qI-$)  and 
det(qI-$)  denotes  the  determinant  of  the  matrix  (qI-$). 
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In  Equation  (2.4) ,  it  is  important  to  recognize  that 

det(qI-<&)  represents  a  polynomial  of  n  order  in  terms  of  the 

forward  operator  q.  Thus,  the  left-hand  side  of  (2.4)  is  a  n 

order  difference  equation  involving  the  output  (or  measurement) 

vectors  z .  Furthermore,  the  elements  of  the  adjoint  matrix  are 

st 

polynomials  of  no  more  than  (n-1)  order  in  q.  Thus,  the  right- 
hand  side  of  (2.4)  is  a  difference  equation  involving  the  input 
vectors  and  w. 

Example:  To  illustrate  the  preceeding  discussion,  suppose 
the  dynamic  model  is  based  on  planar  motion  at  a  nominally  constant 
velocity  except  for  random  forcing  function.  Suppose,  also,  that 
noisy  measurements  are  available  of  the  position  of  the  target. 

The  state  space  model  has  the  form 


r 

> 


J  > 

\ 

\ 

\ 


/: 


v. 

/ 

\ 


or 


x3 (k+1) 

- 1 

t> 

rt 

O 

O 

x2  (k+1) 

0  10  0 

x3  (k+1) 

0  0  1  At 

x^  (k+1) 

m  m 

0  0  0  1 

x(k+l)  •  $x(k)  +  w(k) 


Xj^(k) 

x2(k) 
x3(k) 
x4(k) 

(2.5b) 


w1(k) 

w2(k) 

w3(k) 

wA(k) 


(2.5a) 
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The  measurement  model  is 


or 


z^k)' 

m 

x1(k) 

+ 

v1(k)' 

z2(k) 

x^(k) 

v2(k) 

10  0  0 
0  0  10 


z(k)  *  Hx(k)  +  v(k) 


x^(k) 

+ 

v1(k) 

x2(k) 

v2(k) 

x^(k) 

x4(k) 

(2.6a) 


(2.6b) 


It  follows  that 

det(qI-$) 


(q-1) 


adj(ql-fl>) 


(q-1)  At 
(q-1)3 
0 
0 


0 

0 


0 

0 


(q-1)3  (q-l)2At 
0  (q-1)3 
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<qi-*)‘ 


At 


At 

(q-l)2 


Using  these  results  in  (2. A),  we  have 


(q-1)  z(k) 


(q-1)3  (q-l)2At  0  0 


(q-l)A  /V;L(k)\ 

\v2(k)/ 


wx(k) 
w2(k) 
w3  (k) 
w4(k) 


Thus,  a  fourth-order  difference  equation  for  the  outputs,  a  third- 
order  difference  equation  for  the  plant  inputs,  and  fourth-order 
difference  equation  for  the  measurement  errors  are  obtained. 

From  the  definition  of  the  system,  there  are  actually  two, 

uncoupled,  second-order  systems  being  considered.  This  fact  is 

2 

reflected  by  noting  that  (q-1)  can  be  eliminated  from  the  model 
to  obtain 


(q-1)  S(k)  -  £(k+2)  -  2£(k+l)  +  z/k) 


w^k+1)  -  w3(k)  +Atw2(k) 


w3(k+l)  -w3(k)+AtwA(k) 


+  v(k+2)  -  2v(k+l) + v(k) 
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(2.7) 


> 


> 


» 


► 


» 


» 


» 


z(k+2)  -  2z(k+l)  +  z(k)  -  v(k  +  2)  +  [Bw(k  +  1)  -  2v(k+l)] 


+  [Cw(k)  +  v(k)  ] 

10  0  0 
0  0  10 

0  At  0  0 

0  0  0  At 

Equation  (2.7)  is  a  two-dimensional  description  of  the  behav¬ 
ior  of  the  measurements  as  a  function  of  the  random  plant  inputs  and 
the  measurement  noise.  Again,  it  is  Important  to  note  that  there  is 
no  dependence  upon  the  state  with  the  result  that  no  state  estimation 
is  required. 

Using  the  input-output  model  (2.4),  predictions  of  future  meas¬ 
urements  can  be  obtained.  Since  the  covariance  of  the  random  noise 
is  known,  the  variances  of  the  prediction  errors  can  be  determined. 
There  are  two  major  problems  that  must  be  solved  to  utilize  (2.4)  for 
prediction. 

(1)  White  noise  estimation:  From  the  data,  estimates  of  the 
noise  variables  appearing  on  the  RHS  of  (2.4)  can  be  made. 
These  estimated  values  are  then  used  in  predicting  future 
output  values.  This  problem  is  discussed  in  Section  3. 


where  B  “ 


(2)  Initiating  the  prediction:  Special  considerations  are 
required  to  initiate  the  predictions  and  the  white  noise 
estimetion.  This  problem  is  discussed  in  Section  4. 


► 


3.0  WHITE  NOISE  ESTIMATION 


The  predictive  model  (2.4)  can  be  rewritten  as 

z(k)  +  A-z(k-l)  +  *•*  +  A  z(k-n) 

J.  n — 

-  B.u(k-l)  +  •••  +  B  u(k-n) 
l  n 

+  C.w(k-l)  +  •  •  •  +  C  w(k-n) 
l  n 

+  v(k)  +  A.v(k-l)  +  •••  +  A  v(k-n)  (3.1) 

l  n 

The  matrices  A^ ,  B^,  C^,  |  i  «  1 ,  2,  ***,  n  }  are  known  and  their 

definitions  follow  by  inspection  of  Equation  (2.4).  For  the  example, 
it  is  apparent  from  (2.7)  that  n-2  and 

Ax  -  -21,  A2»I,  CjjB,  C2-C,  B^  ■  B2  “  0 

As  discussed  in  Appendix  A,  the  model  (3.1)  can  be  simplified 
by  modifying  the  noise  sequence  on  the  RHS  [3].  Further,  the  plant 
control  Inputs  ti(i)  shall  be  neglected  for  the  remainder  of  this  dis¬ 
cussion  since  they  are  regarded  as  unknown  but  nonrandom.  The  noise 
sequences  w(i),  v(i)  can  be  replaced  by  a  single  noise  sequence,  say 
n(l),  having  the  same  statistical  properties  and  with  the  property 
that  the  model  is  invertible.  The  procedure  for  determining  the  modi¬ 
fied  model  is  given  in  Appendix  A.  Thus,  the  matrices  D^,  (i “  1,  2, 
•••,  n)  and  the  covariance  matrix  N  are  assumed  to  be  known  and  (3.1) 


becomes 


where  n(i)  represents  a  zero  mean,  white  noise  sequence  with  co- 
variance  matrix  N.  The  model  (3.2)  provides  the  basis  for  the 
remainder  of  our  discussion. 

Consider  the  problem  of  predicting  z(k),  given  measurements 
z(k-l),  z.(k-2),  ••*.  It  is  well  known  [4]  that  the  estimator  that 
minimizes  the  mean-square  error  is  the  conditional  expectation 

z(k/k-l)  £  E[z(k)/z.(k-l),  z(k-2),  •••  ]  (3.3) 

Using  (3.2),  the  predictor  is  given  by 

£(k/k-l)  -  -A^z (k-1)  -  • • •  -  A^z (k-n)  +  n(k/k-l) 

+  D,n(k-l/k-l)  +  **•  +  D  n(k-n/k-l)  (3.4) 

l —  n- 

As  in  (3.3),  the  notational  convention 

n(j/i)  A  E[n(j)|z(i),  z(i-l),  •••  ] 

is  used  in  (3.4).  From  (3.4),  it  is  apparent  that  estimates  of  the 
white  noise  variables  are  required  to  evaluate  the  predictor.  The 
white  noise  estimation  problem  is  considered  in  the  following  para¬ 
graphs.  , 

The  noise  n(i)  is  zero  mean  and  white.  Thus,  one  sees  that 

n(k/k-£)  -  E[n(k)|  z(k-i),  z (k-1-1) ,  •••  ] 

-  E[n(k)] 

-  o , i>i 
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It  remains  to  consider  n(k-i/k-l),  ial.  Given  a  semi-infinite 
span  of  data,  the  white  noise  variables  can  be  determined  precisely 
since  the  noise  model  is  invertible.  To  verify  this  statement,  one 
solves  the  difference  equation  in  (3.2)  for  n(k)  as  a  function  of  the 
semi-infinite  span  of  measurement  values  z.(k),  z^(k-l) ,  ***.  Formally, 
this  can  be  done  by  defining  the  operator 


Dtq'1)  -  I  +  q"^  +  q_2D2  +  •  •  •  +  q"®!^ 


A(q_1)  -  I  +  -  q”2A2  + 


in  terms  of  the  operator  q 


q  ^n(k)  *  n(k-l). 


Then,  Equation  (3.2)  becomes 


A(q  1)£(k)  •  D(q  1)n(k) 


*  +q'“An 


(3.5) 


Since  D(q  A)  is  invertible,  the  solution  of  (3.5)  for  n(k)  is 


n(k)  -  ^Kq’1)]  A(q"*^)_z( 


"  detW-1)  adJ  D<q"1)A(q'1)£(k) 

But  det  D(q  ^)  is  a  polynomial  of  nC^  order  in  the  delay  operator 
q~2.  This  implies  that  the  Inverse,  1/det  D(q-1),  is  an  infinite  series 
in  q~2.  The  Inver tibility  of  D(q  *)  ensures  the  convergence  of  the 
infinite  series.  Thus,  n(k)  has  been  expressed  in  terms  of  the  semi- 
infinite  collection  of  measurements.  Obviously, 

n(k/k)  -  E[n(k)/*(k),  z(k-l),  •••  ] 


ORINCON 


■v  • 


It  follows.  Immediately,  that 


nCk-i/k)  -  n(k-A) ,  l  >  0 

Assuming  the  availability  of  the  data  z,(i) ,  (l»k-l,  k-2, 
•••),  the  predictor  (3.4)  can  be  expressed  as 

£(k/k-l)  -  A1r(k-1) - A^zCk-n) 

+  D.n(k-l)  +  •••  +  D  n(k-n)  (3.6) 

1  n— 

Consider  the  error  in  this  predictor 

7(k/k-l)  -  z(k)  -  z(k/k-l) 

Comparing  (3.6)  and  (3.2),  it  follows  that 

^(k/k-1)  »  n(k)  (3.7) 

Note  that  (3.7)  provides  a  convenient  means  for  computing  n(k)  upon 
receipt  of  the  measurement  z^(k) .  The  covariance  of  the  error  in  the 
estimate  obtains  immediately. 

E[£(k/k-l)£T(k/k-l)]  -  N  ..  , 

Unfortunately,  a  semi- infinite  number  of  measurements  is 
never  available  and  the  noise  sequence  n(k)  must  be  estimated  by  other 
means.  A  reasonable  estimator  n(k/k)  can  be  postulated  directly 
from  (3.2) 

n  n 

n(k/k)  ^  z(k)  +  EA.z(k-i)  -  ED.n(k-i/k-l)  (3.8) 

i-1  1  i-1  1 

This  estimator  requires  knowledge  of  £(k-i/k-i) ,  (i*l,  2,  ••*,  n). 


16 


ORINCON 


t 


> 


I 


» 


{ 


> 


# 


> 


The  estimator  W(k/k)  can  be  wrltcen  as 

D(q_1)n(k/k)  -  A(q"1)*(k)  (3.9) 

Subtracting  (3.9)  from  (3.5)  yields 

»(q“1)In(k)-n(k/k)j  -  0  (3.10) 

Since  D(q  *)  is  invertible  (i. e. ,  stable),  it  follows  that  the  esti¬ 
mation  error  n(k)  -TT(k/k),  must  vanish  as  k  becomes  large  without 
bound.  Thus,  1[(k/k) ,  as  given  by  (3.8),  provides  a  reaaonable  esti¬ 
mator  of  n(k). 

The  covariance  of  the  error  in  the  estimate  n(k/k)  is  easily 
calculated  using  (3.10).  Initial  conditions  for  the  covariance  of 
the  errors  n(i/i) ,  i -  k-1,  ••*,  k-n,  are  required.  For  this  develop¬ 
ment,  it  shall  be  assumed  that  the  errors  in  the  estimates  are  un¬ 
correlated.  Since  n(k)  is  white  noise  and  since  n(k/k)  converges 
to  n(k),  the  errors  must  tend  to  be  uncorrelated.  From  (3.10),  the 
error  ii(k/k)  can  be  written  as 

n(k/k)  -  -D,n(k-l/k-l)  -  •••  -D  n(k-n/k-n) 

—  1  n 

If  the  initial  conditions  are  chosen  to  be  unbiased 


E[n(k-i/k-i)]  -  0,  i-1,  2,  •••,  n 

then  n(k/k)  will  be  unbiased  for  all  k.  Thus,  the  'eovariance  of  the 
error  is 

Etn(k/k)KT(k/k))  -  £  D4E[n(k-i/k-i)KT(k-i/k-i)] 

i-1  1 


P  (k/k) 
n 


L  D.P  (k-i/k-DD.1 
i-1  1  n  1 


Remember  that  P  (k/k)  — ►  0  as  k  — ►  0. 
n 


17 


ORINCON 


1  > 


1 


l  f  » 
f  * 


Consider,  now,  the  error  In  the  predictor  using  7  in  place 
of  n  in  (3.6) 


mm  mm 

s(k)  -  i(k/k-l)  -  -EA.z(k-i)  +  E  D.n(k-i) 

i-1  1  i-0  1 

n  n 

+  L  A.jz  (k-i)  -  L  D  n(k-i/k-i) 


i-1 


i-1 


-  n(k)  +  E  D.  (n(k-i)  -  n(k-i/k-i)  ]  (3.12) 

i-1  1  ” 

From  the  assumptions,  it  is  clear  that  z^(k/k-l)  is  an  unbiased  esti¬ 
mator.  The  covariance  of  the  error  is  given  by 

E  |  (z(k) -£(k/k-l)]  [z(k) -£(k/k-l)]T|>  -  N  +  Pn(k/k)  (3.13) 

The  definition  of  the  predictor  is  now  complete  except  for 
the  consideration  of  initial  conditions  for  the  predictor.  This 
aspect  is  considered  in  Section  4. 


4.0  TRACK  INITIATION 


Suppose  that  the  process  is  initiated  with  the  measurement 
But  the  predictor 

£(k/k-l)  ■  -  j^A  z(k-i)  +  S D .n(k-i/k-i)  (4.1) 

i-1  1  i-1  1 

requires  n  measurements  for  the  one-stage  prediction.  It  also  re¬ 
quires  n  white  noise  estimates. 

In  the  absence  of  the  requisite  data,  the  available  infor¬ 
mation  reduces  to  the  a  priori  statistical  description  of  the  random 
variables.  By  definition, 

E[z(k)]  -  0  for  all  k 
E[n(k)]  ■  0  for  all  k 

The  covariance  of  the  noise  is  known  to  be  N.  The  covariance  of 
the  measurement  js(k)  is  determined  from  the  state  space  model  (2.1)- 
(2.2) 

Z(k)  -  E[z(k)zT(k))  -  HE[x(k)xT(k)]HT  +  R  (4.2) 

where 

E[x(k)xT(k)]  -  $E[x(k-l)xT(k-l)]0T  +  Q 

let 

\  -  E[x(k)xT(k)] 

Note  that  it  satisfies  the  difference  equation 
“k  “  ^“k-l*1  +  Q 


where  an  initial  condition,  say  ,  must  be  specified 
E[x(1)xT(1)]  -  Mx 

Consider  the  prediction  of  z/2)  from  ^(1) .  Then,  (4.1)  re¬ 
duces  to 

*(2|l)  -  -Ajzd)  +  DjSGLll)  (4.3) 

This  represents  an  unbiased  estimate  where  the  a  priori  expected 
values  of  ^(k-i),  n(k-i/k-l) ,  (i-2,  3,  #,,f  4)  have  been  used. 

The  error  in  this  predictor  is 

n  n 

1(2 |l)  -  -Ea.£(2-1)  +  E  D .n(2-i)  +  A.z(l)  -  D.n(l|l) 
i*l  1  i*l  1  1“ 

-  -Ea.z(2-1)  +  E  D.n(2-i)  +  D. [n(l)  -  n(l|l)] 
i-2  1  i-2  i_  1  ~  “ 

The  covariance  is  computed  to  be 

n 

E[i(2|l)IT(2|l))  -  E  A.  Z(2-i)A.T 

i-2  1  1 

n 

♦  Ed.HD.T+  D.P  (l|  1)D.T  (4.4) 

i«2  1  *  1  ®  1 

The  noise  n(l)  au».  be  estimated.  Referring  to ’(3.8) ,  the 
estimate  is 

n(l|l)-z(l)  (4.5) 

with  covariance  matrix 

n  n 

P  (1|1)  -  K[*(1|1)£T(1|1)]-  E  A.Z(l-i)A.T+ED.ND.T  (4.6) 
n  i»l  x  x  i«i  i  i 

This  calculation  completes  the  first  computational  cycle  and  the 
second  measurement  can  be  processed  when  obtained. 
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Suppose  z\l)  has  been  obtained  and  consider  the  prediction 
of  z,(3).  The  noise  sample  n(2)  is  estimated,  using  (3.7)  as 


n(2 | 2)  -  z (2)  -  1(2 |l) 


with  covariance  matrix 


P  (2 1 2)  -  £  A.Z(2-i)A.T+D.P  (IIDdJ+Ed.ND. T 
n  i«2  A  n  A  i-2  1  x 


The  predictor  is  given  by 


£(3 1 2)  -  -A12(2)-A22(l)+D1n(2|2)+D2n(l|l) 


The  predictor  has  the  covariance  matrix 


(4.7) 


(4.8) 


(4.9) 


E[z;(3|2)!T(3|2)]  -  E  A  Z(3-i)A?  +  £d.P  (3-i/3-i)D.1 

i-3  1  x  i-1  1  n  i 


+  Z-D.ND, 


(4.10) 


Equations  (4.7)  -  (4.10)  define  the  second  computational  cycle. 

The  algorithm  described  for  k»l  and  k ■  2  generalizes  read¬ 
ily  for  k  £  n.  The  general  procedure  is  stated  below.  For 
k-1,  2,  '•*,  n 


n(k/k)  ■  z(k)  -  z(k/k-l) 


M  U 

P  (k/k)  -  E  A.Z(k-i)A,T  +  ED,ND,T 
n  i-k  1  1  i-k  1  1 

k-1 

+  E  D.P  (k-i/k-i)D .T 

i»l  i  n  i 

k-1  k-1 

z (k/k-1)  -  -E  A.z(k-i)  +  E  D.n(k-i/k-i) 

i-1  i_  i-1  i_ 


(4.11) 


(4.12) 


(4.13) 
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where  2(1 |0)  ■  0. 


The  system  (4.11)-  (4.14)  applies  for  k>n  by  following  the 

n  k-1  n 

convention  that  summations  ^  ( • )  »  0  for  k  >  n  and  ^  ( * )  *  £  ( * )  • 

i**k  i"*l  i“l 

Thus,  the  algorithm  is  conpletely  defined  by  (4.11)  -  (4.14)  for  all 


k.  The  covariance  matrix  Z(k)  is  given  by  (4.2). 


5.0 


MULTITARGET  TRACKING  ALGORITHM 


The  prediction  algorithm  (4.11)-  (4.14)  provides  the  basis 
for  data  association  for  the  multitarget  tracking  algorithm.  In 
the  following  discussion,  the  data  association  algorithm  is  pre¬ 
sented.  The  predictions  utilized  in  the  algorithm  are  assumed  to 
be  obtained  from  the  algorithm  (4.11)  -  (4.14) .  The  simplicity  of 
this  algorithm  provides  a  significant  computational  advantage  rela¬ 
tive  to  existing  algorithms  which  utilize  the  Kalman  filter  algo¬ 
rithm  to  accomplish  the  predictions.  No  filter  calculations  are 
required  here  for  track  prediction.  The  track  predictions  are 
carried  out  in  terms  of  output  variables  which  are  generally  of 
lesser  dimension  than  the  state  dimension.  Having  established 
feasible  tracks  using  the  input-output  model,  the  Kalman  filter 
can  be  used  for  state  estimation  for  only  the  feasible  tracks. 

Only  the  data  association  problem  is  addressed  in  the  following 
paragraphs. 

Consider  a  measurement  time  t^  and  suppose  measurements 

*••,  Zj^t^)}  are  obtained.  At  t^^*  suppose  that 
there  are  ^  track  files.  A  track  file  is  regarded  as  a  collec¬ 
tion  of  data  and  white  noise  estimates  that  are  used  to  predict  the 
next  measurement  sample  for  the  track.  These  track  files  shall  be 
denoted  as  f^_^^(z^  \  n^  S •  (i*l,  2,  •••,  i^^) .  EaCh  track  file 
can  be  used  to  generate  a  measurement  prediction  z^(k/k-l)  with  asso¬ 
ciated  covariance 

P  ±(k/k-l)  -  E[z\(k/k-l)?.T(k/k-l)] 
z  — 1  — X 

The  corresponding  error  ellipsoid  can  be  used  to  define  a  data  win¬ 
dow.  Measurements  z, (t.  )  contained  within  the  window  are  associated 

— 1  k 

with  the  track  file.  For  each  track  prediction,  one  of  three  occur¬ 
rences  is  possible. 
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(1)  No  measurements  z^(t^)  are  contained  In  the  error  ellip¬ 
soid  of  the  j6*1  track  file. 

In  this  case,  one  can  continue  the  track  file  by  predicting 
to  the  next  measurement  time  and  repeating  the  test.  Although  only 
one-stage  prediction  has  been  discussed,  the  algorithm  (4.11)  -  (4.14) 
is  adapted  in  a  simple  manner  to  accomplish  N-stage  predictions. 

Alternatively,  the  absence  of  feasible  measurements  may  indi¬ 
cate  the  infeasibility  of  the  track  file  with  the  result  that  the 
file  is  dropped  from  memory.  The  logic  for  establishing  conditions 
for  eliminating  specific  track  files  is  dependent  upon  statistical 
thresholds  that  are  defined  by  the  system  designer.  The  objective 
of  data  association  is  to  reduce  the  number  of  potential  target 
tracks.  Thus,  it  is  imperative  that  track  files  are  eliminated 
whenever  reasonable. 


(2)  One  measurement  (t^)  is  contained  in  the  error  ellip¬ 
soid  of  the  track  file. 

This  circumstance  provides  clear  indication  of  the  need  for 
continuing  the  track  file.  The  measurement  is  added  to  the  file 
and  the  track  file  is  updated.  Note  that,  generally,  the  data 
z^(tk_n_j)  is  not  needed  for  prediction.  It  may  be  necessary  to 
retain  this  measurement  sample  for  input  to  the  state  estimator  in 
the  event  that  this  track  file  is  identified  as  being  feasible. 

It  can  happen  that  the  measurement  lies  in  the  data 

window  for  more  than  one  track  file.  From  physical  considerations, 
a  measurement  must  relate  to  a  single  target.  However,  additional 
data  may  be  required  before  a  unique  assignment  of  _z^(t.  )  can  be 
achieved.  Alternatively,  the  measurement  may  be  arbitrarily  assign¬ 
ed  at  t^  to  a  specific  track  based  on  other  considerations.  For  ex¬ 
ample,  it  may  be  assigned  to  the  track  yielding  the  smallest  predic¬ 
tion  error  [z^ (t^)  - %.  (t^/t^^) ]  •  Alternatively,  it  may  be  assigned 
to  the  track  having  the  smallest  error  ellipsoid  volume. 


(3)  More  than  one  measurement  is  contained  in  the  data  win¬ 
dow  for  the  track  file. 

At  least  two  alternatives  present  themselves  in  this  situation. 
First,  the  measurement  yielding  the  smallest  prediction  error  (e.g. , 

T 

smallest  in  the  sense  that  the  squared-residual  -Jy 

[z^(t^)  -  Sj (t^/t^-j) 1  is  least  for  all  z^(t^)  contained  in  the  data 
window)  may  be  added  to  the  track  file  and  all  other  measurements  re¬ 
jected.  Alternatively,  a  splitting  of  track  files  may  be  introduced._ 
The  track  file  f ( • , • )  may  be  used  to  generate  additional  track 
files  which  differ  only  through  the  measurements  z^(t^)  which  are 
contained  in  the  data  window. 

Some  measurements  z.  (t. )  might  not  be  contained  in  the  data 
window  for  any  track  file  or  might  not  be  included  in  a  track  file 
when  track  splitting  is  not  permitted.  It  is  important  that  every 
measurement  is  considered  as  emanating  from  a  potential  track.  Meas¬ 
urements  not  added  to  a  track  file  should  be  used  to  initiate  a  new 
track  file.  As  discussed  in  Section  4,  the  prediction  for  tracks  con¬ 
taining  fewer  measurements  than  the  order  of  the  model  uses  a  priori 
statistical  information  about  the  measurement  and  input  noise  covar¬ 
iances.  Care  must  be  taken  in  defining  these  covariances  in  order 
to  prevent  the  data  windows  for  subsequent  predictions  from  being 
unreasonably  large. 


From  the  preceding  discussion ,  one  notes  that  at  t^  there 
will  be  at  least  M,  track  files.  Generally,  there  will  be  many  more 


track  files. v  As  this  could  be  an  enormous  number,  it  is  im¬ 
perative  that  infeasible  tracks  be  eliminated  as  early  as  possible. 
This  elimination  occurs  when  no  measurements  are  contained  in  a  track 
data  window  for  one  or  more  sampling  times. 

A  feasible  track  is  Identified  when  the  number  of  measurements 
in  a  track  file  becomes  sufficiently  large  (i.e.9  k  >  n) .  When  a 


track  is  identified,  the  file  is  retained  as  long  as  relevant  data 
is  obtained.  But,  also,  the  data  can  be  provided  to  a  state  estima¬ 
tor.  In  addition  to  obtaining  estimates  of  the  state  of  the  track, 
the  feasibility  of  the  track  can  be  confirmed  or  refined  by  analysis 
of  the  resulting  innovations  sequence.  This  possibility  shall  only 
be  noted  here  but  shall  not  be  pursued  further. 

A  general  algorithm  has  been  defined  for  accomplishing  data 
association  and  feasible  track  identification  for  the  multitarget 
tracking  problem.  For  the  class  of  linear,  time-invariant  systems 
that  have  been  considered,  it  appears  that  this  procedure  may  be 
considerably  more  computationally  efficient  than  procedures  reviewed 
in  Reference  2  which  are  based  upon  the  use  of  the  Kalman  filter. 

The  utility  of  the  procedure  depends,  as  for  the  KF-based  methods, 
upon  the  definition  of  data  windows  that  permit  the  elimination  of 
infeasible  tracks  in  a  statistically  satisfactory  manner.  The  need 
for  state  estimation  to  achieve  data  association  has  been  eliminated 
and  a  simpler  prediction  using  equivalent  input-output  models  has 
been  proposed. 

Because  of  the  inherent  complexity  of  the  multitarget  track¬ 
ing  problem,  further  analysis  of  this  algorithm  must  be  accomplished 
using  numerical  experiments.  The  detailed  computations  required  for 
this  type  of  analysis  could  not  be  performed  due  to  a  lack  of  funds 
remaining  in  this  contract.  The  approach  that  has  beert  defined 
appears  to  be  extremely  promising.  Additional  contractual  support 
to  Investigate  this  approach  in  greater  depth  seems  warranted. 
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APPENDIX  A:  SIMPLIFYING  THE  NOISE  MODEL 

The  noise  terns  in  the  RHS  of  (3.1)  involves  two  white  noise 
sequences.  To  simplify  the  development  of  the  predictor,  it  is  de¬ 
sirable  to  introduce  a  separate  noise  sequence,  say  n(i),  that  yields 
the  same  noise  characteristics  as  w(i)  and  v(l).  That  is,  we  want 
to  determine  coefficient  matrices  such  that 

n(k)  +  D^n(k-l)  +  •  •  •  +  D^r|(k-n)  -  v(k)  +  A^k-l) 


+  •••  +  A  v(k-n)  +  C.w(k-l)  +  •••  +  C  w(k-n) 
rr-  X“  nr- 


c(k)  =  cx(k)  +  ;2(k) 


(A.l) 


-  eJ [  £  A^k-ljj  [  £  A^lk-l-jtj  | 


1 

E[?2(k)i2T0t-j)]  -  eI[£^H[£  C^Ck-i-j  )  j 


E  C.QCj 


i«3+l 


ixwi-j 


For  the  new  sequence  £(k),  it  follows,  since  w  and  v  are  in¬ 
dependent  ,  that 

E[£(k)£T (k-J  )  ]  -  E tij (k)^1  (k-n)  ]  +  E(£2  (k)£2T (k-n)  ] 


-  £aX  a  +  Ec.QC* 

i-j  1  1  J  i-j+1  1  1  J 


0  for  all  j  >  n 


(A.  2) 


But  we  want 
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i-j 


DjND 


T 

i-j 


-  0  for  j  >  n  (A.  3) 

Equating  coefficients  in  (A. 2)  and  (A. 3) ,  one  obtains  (n  + 1)  matrix 
unknowns  D^,  D2,  Dq,  N. 

D  N  -  A  R 
n  n 


D  ,N  +  D  ND 
n-1  n 


X 


A  ,R 
n-1 


E  D  ND  T  -  E  A  RA  T  + 
i-0  1  1  i-0  1  1 


V*  x 
E  C.QC. 
i»l  i  1 


(A.  4) 


All  of  the  terms  on  the  RHS  of  (A. 4)  are  known.  It  is  possi¬ 
ble  to  solve  for  the  D^,  i»l,  2,  •••,  n  and  for  N.  More  than  one 
solution  is  possible,  generally.  Then,  the  solution  is  chosen  to  in¬ 
sure  that  the  system  is  invertible  in  the  sense  that  the  zeros  of 
the  determinant  of 

D(q)  -  qnl  +  +  •••  +  Dq 

are  all  contained  within  the  unit  circle.  With  the  satisfaction  of 
this  condition,  the  difference  equation 

D(q)n(k)  -  0 

is  stable  and  n(k)— ♦O  as  k  — ►  00 . 
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EXAMPLE 


As  derived  earlier,  consider  the  noise  model  given  by 
ijOO  +  £2(k)  -  v(k)  -  2v(k-l)  +  v(k-2)  +  Bv(k-l)  +  Cv(k-2) 

where  n  -  2  and 


Consider  (A. 4).  This  system  is  given  below. 
d2n  -  a2r  -  R  -  I 
DjN  +  D2ND1  -  AjR  +  C2QC1T  -  -2o2I 

N  +  D^NDj1  +  D2OT>2T  -  R  +  Aj^RAj^1  +  A2RA2T  +  CjQCj1  +  CjQCj1 
-02I  +402I  +02I  +Ou2l  +  At2oM2; 

"  H  +  °U  +  At2°222  ]  1 

Substituting  D2N  into  the  second  and  third  equations,  we 

obtain 

D,N  +  a2D-  -  2o2  I 
1  v  1  v 

T  2  T  T  2  2  2  21 

N  +  DjND^  +  0‘  -  [6oy  +On  +  At  022  J  I 
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N 2  (  <J4  1+2 O2  N  +  N2)  +4o2N2+a4  (o4  I  +  2o4  N  +  N2) 
\  v  V  /  V  v\v  V  / 

2  2  \ 

I  +  2ov  N  +  N  N 

4  23  22  22  42  8  8'** 

N  +2o  N  +0,  N  +4a  N^  +  a  N  +2o  N  +  O  I 
v  4  v  v  v  v 


3  2  2  4 

-  aN  +  2ao  N  +  ao  N 
v  v 


ORINCON 


This  equation  can  be  satisfied  by  assuming  N  to  have  the 

form 

N  -  nl 

Then,  a  scalar,  quart ic  equation  must  be  solved  to  complete  the 
definition  of  the  modified  system. 

"*+  K“)'>3+  M^2+  K-*)  v+°v-° 

The  value  of  N  is  used  to  determine  and 
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A  DISCUSSION  OF  THE  OCEAN 
TACTICAL  TARGETING  PROBLEM 

1.0  INTRODUCTION 

Any  realistic  approach  to  the  ocean  targeting  problem  requires 
the  interaction  of  many  diverse  sensor  and  computational  systems; 
analytical  and  computational  techniques;  and  many  political  empires. 
The  fusion  of  these  elements  is  fundamental  to  any  solution  to  the 
OT^  problem.  This  white  paper  is  an  attempt  to  describe  the  problem 
from  a  systems  scientific  point  of  view  and  to  highlight  problem  areas. 

Consider  a  large  ocean  surveillance  area  with  well  defined 
boundaries.  This  surveillance  area  contains  surface  and  subsurface 
ships  as  well  as  aircraft  and  other  signal  sources.  These  may  be 
threats,  potential  threats,  friendly s  or  neutrals.  Each  of  these 
categories  could  be  broken  down  into  finer  divisions — such  as 
combatant,  merchant,  fishing — ending  in  hull  number  or  specific 
identification  of  all  sources. 

Resources  including  sensors,  computational  resources, 
communication  resources  and  weapon  resources  are  available  to  help 
monitor  and  control  the  ocean  surveillance  area.  There  are  numerous 
sensors  and  sensor  systems  that  gather  data  containing  information 
about  targets  in  the  surveillance  area  or  about  passages  into  or  out 
of  the  area.  There  are  communication  links  that  allow  transmission 
of  raw  and  processed  sensor  data;  information  about  potential 
target  tracks  and  identification;  and  information  used  to  control  or 
allocate  sensors,  computational,  personnel  and  weapon  resources. 

There  are  computers  used  to  perform  the  computational 
tasks  of  signal  processing,  data  association,  detection,  tracking, 
classification  and  resource  allocation.  There  are  weapons  to  be  allo¬ 
cated  to  threat  targets  and,  last  but  not  least,  personnel  to  be 
allocated  to  various  aspects  of  the  problem  requiring  human  interaction 
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Thus  the  ocean  targeting  problem  is  the  problem  of  utilizing 
these  available  resources  to  first  obtain  an  accurate  picture  of  the 
surveillance  area  and  secondly  to  allow  adequate  weapon  allocation 
to  be  made.  This  can  be  visualized  as  a  classical  control  system 
with  two  feedback  loops  as  shown  in  Figure  1.  The  feed- forward 
loop  consists  of  the  surveillance  area  and  all  that  it  contains:  the 
sensors,  sensor  system  processing  and  fusion  processing  required  to 
do  classical  ocean  surveillance. 

The  first  or  inner  loop  is  the  sensor  control  loop  which 
allows  control  of  the  sensor  resources.  This  allows  the  focusing 
of  the  sensor  resources  on  a  particular  region  of  interest  or  a  par¬ 
ticular  target  of  interest.  Thus  sensors  can  be  focused  on  a  particu¬ 
lar  geographical  or  frequency  region  or  on  events  of  interest  because 
of  a  priori  information,  or  information  obtained  from  completely  differ¬ 
ent  sensor  systems.  The  utilization  of  this  inner  loop  in  an  intelligent 
manner  is  vital  to  the  success  of  an  ocean  targeting  system  but  it 
also  implies  great  problems  of  a  theoretical,  computational  and  political 
nature.  This  is  one  of  the  major  unsolved  problems  in  development 
of  a  real  ocean  targeting  system. 

The  second  feedback  in  Figure  1  is  the  weapon  allocation /fire 
control  loop.  This  is  the  major  role  of  command,  control  and  commu¬ 
nication  systems.  It  should  be  clear  that  any  design  of  a  weapon 
allocation  system  will  have  great  impact  on  the  sensor  allocation 
loop  and  the  processing  done  in  the  feed-forward  ocean  surveillance 
path.  In  order  to  reduce  the  complexity  of  this  discussion,  we  will 
first  focus  on  the  feed-forward  path  of  the  ocean  tactical  targeting 
problem  without  considering  the  feedback  paths.  The  effects  and 
problems  introduced  by  the  inner  and  outer  feedback  loops  will  be 
considered  later. 

The  difficulties  inherent  in  the  feed-forward  path  or  in  the 
surveillance  problem  can  be  broken  down  into  several  areas.  These 
include  the  standard  multi-target,  multi-sensor  detection,  tracking. 
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and  classification  problem  in  clutter;  the  data  communication  problem; 
the  distributed  processing  problem;  the  data  base  problem;  the  man/ 
machine  interface  problem;  and  the  political  problem. 

A  somewhat  more  complex  view  of  the  OT^  problem  is  shown 

in  Figure  2.  Here  many  possible  levels  of  sensor  system  signal 

processing  are  shown.  Also,  both  the  inner  and  outer  feedback 

possibilities  from  Figure  1  are  indicated.  Note  that  the  inner 

sensor  control  feedback  path  is  really  a  distributed  feedback  loop. 

Thus  at  least  the  possibility  exists  for  sensor  level  control  at  each 

sensor,  each  sensor  system  level  and  at  the  overall  fusion  center. 

3 

Sensor  control  could  conceptually  be  exercised  from  the  C  function 
in  order  to  improve  fire  control  or  targeting  functions. 

There  are  many  difficulties  inherent  in  the  feed-forward  path 
or  in  the  surveillance  problem.  Some  of  these,  which  ORINCON 
feels  are  key  research  areas,  are  discussed  in  the  sections  to 
follow . 

Section  2  is  a  discussion  of  the  overall  ocean  tactical  targeting 
problem  with  particular  consideration  of  various  data  partitions. 

In  Section  3  the  problems  and  potential  solutions  to  data 
communication  will  be  addressed.  Two  approaches  appear  viable; 
namely,  data  compression  and  sensor  level  preprocessing. 

Section  4  presents  the  key  features  pertaining  to  distributed 
information  processing  as  it  relates  to  ocean  tactical  targeting. 

It  discusses  a  host  of  scenarios  that  could  be  implemented. 

A  number  of  problems  in  data  base  management  will  be 
illustrated  in  Section  5.  It  then  goes  on  to  describe  certain  tech¬ 
niques  that  could  be  applied  to  attack  these  problems. 

In  Section  6,  a  general  discussion  of  the  measurement  fusion 
problem  is  presented.  Following  this,  in  Section  7,  the  multi-sensor/ 
multi-target  tracking,  detection,  and  classification  problem  is  treated. 
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One  of  a  number  of  possible  approaches  is  discussed  in  some  detail 
to  illustrate  the  complexity  of  the  measurement  fusion  problem. 

Finally,  in  Sections  8  and  9*  the  issues  of  force  and  weapon 
allocation  and  final  targeting  will  be  addressed. 
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2.0  THE  OCEAN  SURVEILLANCE  PROBLEM 


The  ocean  surveillance  problem  as  defined  in  the  last  section 
is  to  take  and  process  data  from  all  available  sensors  which  contain 
information  about  the  ocean  surveillance  area  or  passages  through  its 
boundary,  and  produce  the  best  possible  picture  of  the  surveillance 
area.  This  requires  the  detection,  tracking  and  classification  of  all 
possible  targets  and  the  allocation  of  the  computational,  communica¬ 
tion  and  personnel  resources  available.  This  is  summarized  in  an 
extremely  simple  manner  in  Figure  2,  which  shows  the  possible  sensors 
arranged  vertically.  Any  spatial  factors  such  as  beams,  time  delay 
and  Doppler  resolution  cells,  range  and  range  rate  cells  or  positional 
cells  (of  a  particular  sensor)  are  indicated  by  the  horizontal  "spatial 
sector"  direction.  The  time  flow  of  data  for  a  given  sensor  and 
spatial  sector  (perhaps  a  beam)  is  indicated  by  the  axes  into  the 
paper.  Thus  this  figure  simply  indicates  that  the  overall  "open  loop" 
surveillance  problem  consists  of  taking  all  the  data  from  a  suite  of 
sensor  systems  over  a  time  interval  and  drawing  an  accurate  picture 
of  the  surveillance  area.  This  accurate  picture  of  the  surveillance 
area  must  show  the  time  evolution  of  the  states  of  all  potential  targets. 

It  requires  that  the  functions  of  target  detection,  target  tracking  and 
target  identification  be  done,  as  well  as  possible,  from  the  sensor  data 
available.  The  target  state  will  include  position  and  velocity  information 
as  well  as  any  other  information  required  to  specify  the  future  state  of 
the  target  based  on  the  current  state  and  assumed  dynamics.  This 
could  include  frequency  information,  local  sound  speed  information,  and 
tums-per-knot  information,  as  well  as  parameters  of  the  system 
dynamics. 

In  actual  practice  it  is  impossible  to  take  all  of  the  information 
from  the  sensors  and  process  it  simultaneously  to  obtain  an  overall 
picture  of  the  surveillance  area  or— as  in  Figure  2— an  "all  sensors 
master  track  file."  However,  it  is  clear  that  if  this  were  possible, 
the  maximum  information  could  be  extracted  from  the  raw  sensor  informa¬ 
tion,  and  the  best  possible  system  performance  could  be  obtained. 
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In  any  realistic  attempt  to  solve  this  problem*  intermediate 
processing  steps  must  be  performed.  This  is  required  to  break  the 
problem  down  into  easily  handled  segments  and  is  intended  to  reduce 
the  overall  amount  of  processing  required*  to  reduce  the  amount  of 
data  transmitted,  and  to  allow  the  staff  working  at  a  given  sensor 
level  to  monitor  the  progress  of  the  system.  This  will  generally 
lead  to  performing  operations  at  a  local  level  when  possible. 

Thus,  a  realistic  approach  to  the  problem  might  be  as  indicated 
in  Figure  3.  Here  data  from  each  sensor  system  is  processed 
repeatedly  to  form  a  master  track  file  for  that  sensor  system.  This 
overall  track  information  could  then  be  merged  in  a  fusion  facility  at 
the  track  level  to  form  an  all’- sensor  master  track  file.  It  should  be 
noted  that  information  could  be  lost  in  such  a  system.  If  two  pieces 
of  information  from  two  different  sensor  systems  were  sufficient  to 
form  a  good  track  when  put  together  but  not  sufficient  to  give  tracks 
in  the  separate  sensor  track  files,  then  the  track  would  be  lost.  None¬ 
theless,  some  loss  of  information  will  probably  have  to  be  tolerated 
in  order  to  allow  the  global  problem  to  be  handled  in  even  a  sub- 
optimal  manner. 

In  actual  practice  the  current  system  looks  more  like  Figure  4, 
where  tracking  is  done  on  the  output  of  each  individual  sensor  (per¬ 
haps  in  sensor  coordinates) ,  and  then  all  of  the  tracks  in  sensor 
coordinates  were  put  together  to  form  a  sensor  system  master  track 
file.  The  last  step  of  combining  the  tracks  from  numerous  sensor 

systems  is  done  in  a  very  crude  manual  fashion. 

% 

Note  that  the  possibility  for  all  of  these  levels  of  processing 
was  contained  in  the  more  complex  view  of  the  OT^  problem  shown 
in  Figure  2.  There  many  of  the  possible  levels  of  sensor  system 
processing  are  shown  as  in  Figure  4. 
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3.0  DATA  COMMUNICATION  PROBLEMS 

The  data  communications  problem  is  of  paramount  importance 
in  a  large,  distributed  surveillance  network.  The  raw  information 
rates  produced  by  sensor  systems  in  the  network  will  easily  overwhelm 
the  relatively  low  capacity  transmission  links  that  connect' the  network 
nodes.  Thus  the  network  structure,  and  particularly  the  distribution 
of  processing  resources,  will  be  fundamentally  constrained  by  the 
data  communication  resources  available.  The  two  principal  qualities 
of  these  resources  are  capacity  and  vulnerability. 

Network  link  capacity  is  purchased  at  a  premium  price  in  ocean 
surveillance.  Dedicated,  long-haul  channel  bandwidth  is  very  expen¬ 
sive,  whether  it  be  a  land-based  or  satellite  link.  There  are  two 
approaches  to  minimizing  the  capacity  requirements:  data  compression 
and  sensor-level  preprocessing.  The  data  compression  alternative 
involves  spatial  and/or  temporal  encoding  of  the  outputs  of  a  given 
sensor  system,  so  that  they  may  be  transmitted  at  a  substantially 
reduced  rate,  and  reconstructed  with  acceptable  fidelity  at  the 
receiving  node.  The  advantage  of  this  method  is  that  full  data  pro¬ 
cessing  flexibility  is  preserved  at  the  receiving  node,  and  therefore 
the  maximum  amount  of  information  may  be  extracted  from  the  sensor 
system,  particularly  in  the  case  of  multi- sensor  processing. 

The  sensor-level  preprocessing  alternative  involves  a  more 
extensive,  and  basically  irreversible,  reduction  in  the  information 
rate  of  the  sensor  through  the  extraction  of  certain  parameters  or 
features  of  the  output  (e.g.,  detections,  classifications,  track  infor¬ 
mation).  Only  this  information  is  transmitted  over  the  link,  at 
extremely  low  rates  relative  to  the  raw  data  rate.  The  latter  advan¬ 
tage  is  offset  by  the  reduced  flexibility  of  the  received  data,  i.e., 
further  processing  options  are  restricted  by  the  much-reduced  amount 
of  information  available.  Also,  the  individual  sensor  systems  will 
tend  to  be  more  complex  due  to  their  processing  capabilities,  and 
therefore  more  costly.  It  is  likely  that  both  of  the  alternatives  that 
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have  been  discussed  would  be  incorporated  into  various  portions  of 
the  ocean  surveillance  network. 

Network  vulnerability  includes  considerations  of  data  encryption, 
link  jamming  resistance,  and  error  correction  capabilities.  Data 
encryption  of  both  sensor  outputs,  as  well  as  command  and  control 
information,  is  important  to  the  security  and  effectiveness  of  the 
surveillance  network.  Depending  on  the  political  situation  in  some 
cases,  covert  transmission  using  spread-spectrum  techniques  may 
be  required.  These  techniques  may  also  be  necessary  to  ensure  reliable 
transmission  from  deployed  sensor  systems  that  are  susceptible  to  jam¬ 
ming  tactics  by  the  enemy.  Finally,  the  noisy  characteristics  of 
certain  links,  such  as  telephone  lines,  may  dictate  the  use  of  error 
correction  and  control  methods,  so  that  data  "drop-outs"  can  be 
avoided. 

In  addition  to  considering  the  technology'  available  for  network 

links,  several  issues  regarding  the  data  communications  network  as 

a  whole  must  be  considered.  The  design  of  a  sensor  and  control  data 

network  is  very  sensitive  to  the  particular  solution  used  for  the 
2 

general  OT  problem.  For  instance,  if  the  solution  approach  to  dealing 
with  sensor  data  involves  primarily  gathering  data  into  a  centralized 
fusion  site  for  processing,  then  the  network  topology  would  be  a 
centralized  topology — a  larger  version  of  what  is  used  for  gathering 
data  for  the  ARC. 

However,  if  issues  of  survivability  enter  in — as  they  should — 
then  other  network  topologies  would  be  more  appropriate.  If  a  satellite 
network  was  to  be  used  for  the  network  backbone,  then  a  more 
desirable  topology  class  would  be  the  star  polygon.  These  topologies 
have  several  advantages.  Most  importantly,  enemy  action  cannot  take 
out  a  central  node — thus  destroying  the  network.  Furthermore,  this 
approach  is  convenient  when  a  worldwide  distribution  of  sources  and 
sinks  is  required.  An  economic  aspect  is  that  many  satellites  with 
identical  hardware  may  be  used — and  if  the  satellites  are  mission- 
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reconfigurable  by  uplinked  software,  then  the  network  may  be  quickly 
(even  automatically)  adapted  to  new  satellites  or  satellites  taken  out 
of  service. 
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Any  network  that  supports  a  global  distribution  of  sources  and 
sinks  would  provide  the  cross-linking  of  sensor  data  and  control 
messages  necessary  to  fully  utilize  concepts  of  distributed  processing 
and  distributed  data  bases.  However,  any  network  of  that  scope 
would  clearly  have  to  support  multiple  missions  and  a  wide  range  of 
message  traffic  with  different  bandwidths,  burst  characteristics, 
and  destination  distribution  (one  sink  versus  many  sinks) .  Several 
of  the  more  common  network  management  schemes  are  dedicated  circuit, 
circuit  switching,  packet  switching,  and  message  switching.  Each  of 
these  methods  has  advantages  for  certain  classes  of  message  traffic. 
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There  are  several  areas  of  investigation  appropriate  to  the 

network  issue.  One  area  is  determining  the  utility  of  various  data 

2 

network  topologies  applied  specifically  to  OT  .  Associated  with  this 
area  is  the  question  of  what  mix  of  management  technologies  (packet 
switching,  circuit  switching,  etc.)  is  appropriate  for  OT^  communica¬ 
tion  requirements.  If  multiple  users  are  involved,  how  do  many  users 
access  a  common  system  given  the  need  for  security  and  anti-jamming 
methods?  Since  different  communications  technologies  are  involved 
(satellite,  land  line,  microwave,  etc.),  what  effects  will  data  lateness 
over  various  links  of  the  network  have  on  a  specific  user  such  as  a 
fusion  center?  This  is  especially  important  since,  if  adaptive  routing 
of  messages  is  used  to  accommodate  traffic  jams  or  link  failures  (like 
in  the  ARPANET) ,  data  might  have  random  arrival  times  as  opposed 
to  predictable  arrivals.  Finally,  the  performance  of  the  network  under 
attack  and  crisis  scenarios  should  be  studied.  A  mix  of  analytic  and 
simulative  techniques  will  probably  be  necessary  to  identify  network 
problems  due  to  missing  or  degraded  components. 
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4.0  DISTRIBUTED  INFORMATION  PROCESSING 


Any  solution  to  the  Ocean  Tactical  Targeting  problem  will 
require  digital  and  analog  processing  techniques  at  all  levels. 

Examples  of  digital  processing  tasks  at  the  lower  levels  are  prefilter¬ 
ing  and  calibration  of  data  at  the  sensor  level,  data  compression  or 
encoding  of  data,  and  sensor  level  data  analysis.  At  the  higher  levels, 
some  tasks  are  management  of  data  bases  at  the  tactical  and  fleet 
level,  sensor  data  integration  or  fusion,  target  detection  and  tracking, 
threat  type /target  association,  threat  prioritization  at  the  tactical 
level  (e.g.,  incoming  missiles,  approaching  aircraft),  and  fleet  level 
force /threat  evaluation  and  assignment  (e.g.,  assignment  of  specific 
ASW  components  to  an  identified  submarine  threat) .  One  possible 
approach  to  the  system  design  is  to  ship  data  from  all  sensor  sites 
to  a  central  evaluation  center  where  all  of  the  processing  mentioned 
above  is  carried  out.  Fleet  components  could  be  commanded  based  on 
evaluation  of  the  input  data.  There  are  several  reasons  why  this 
approach  is  impractical  in  general. 

First  and  foremost  is  that  the  data  rate  of  some  sensor  systems 
is  much  too  high  to  be  rooted  from  a  sensor  site  to  an  evaluation 
center.  Some  processing  is  necessary  to  select  a  data  subset,  encode 
or  compress  the  data  to  reduce  the  required  transmission  bandwidth, 
or  to  perform  pre-analysis  of  the  data  to  extract  a  desired  component 
of  the  information  available  in  the  data — again  a  bandwidth  reducing 
function.  At  least  at  the  lowest  level,  then,  certain  digital  and 
analog  processing  tasks  must  be  'distributed'  because  of  communication 
bandwidth  limitations. 

Another  reason  for  distribution  of  processing  functions  is  sur¬ 
vivability  and  flexibility.  As  with  the  data  communications  problem, 
a  level  of  redundancy  is  necessary  so  that  if  part  of  the  system  fails, 
the  remaining  resources  will  (perhaps  at  some  degraded  level)  still 
perform  the  ocean  targeting  function.  An  issue  which  is  closely 
related  to  system  survivability  is  that  of  flexibility.  Since  many 
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elements  of  an  ocean  targeting  system  would  be  mobile  (weapons 
systems  and  sensors  in  both  airborne  and  ocean-based  platforms) ,  a 
fleet  commander  might  wish  to  allocate  a  larger  portion  of  all  resources 
to  a  particular  "hot  spot."  Physically  placing  the  computational 
resources  in  or  near  the  activity  site  could  yield  benefits  such  as 
improved  coverage,  finer  targeting  rsolution  and  shortened  system 
response  time.  Thus  a  commander  could  'zoom  in'  to  a  particular  area 
of  interest  for  special  action. 

A  third  reason  for  distributing  computational  resources  is 
the  concept  of  task  sharing.  Some  system  tasks  may  be  designed  so 
that  they  can  be  processed  at  several  different  processing  nodes. 

If  processing  nodes  in  one  part  of  the  system  become  overloaded, 
tasks  could  be  shifted  to  make  use  of  underutilized  processing 
resources  at  other  nodes.  For  example,  destroyers  in  support  of 
a  carrier  could  process  radar  data  locally  while  linking  sonar  data 
to  the  carrier  where,  presumably,  a  larger  processing  resource 
would  reside.  There  are,  of  course,  problems  with  doing  this  in 
general.  The  first  problem  is  that  some  processing  resources  in 
the  system  have  unique  structures  and  have  permanent  specializations 
(e.g.,  beam  forming  hardware  for  acoustic  data  sources — this 
equipment  can  be  ’retasked'  only  in  the  sense  of  sensor  allocation 
to  specific  geographic  targets) . 

A  second  problem  associated  with  task  sharing  is  the  volume 
of  intemode  communication  traffic  required  to  completely  shift  the 
processing  may  be  impractical.  A  tactical  data  base  could  not  be 
transferred  from  one  ship  to  another  with  surplus  processing  capa¬ 
bility  because  the  lateness  associated  with  the  arrival  of  the  complete 
data  base  might  be  comparable  to  the  flight  time  bf  an  incoming 
missile. 

However,  there  are  many  areas  in  which  this  could  be  a  valuable 
concept.  One  (admittedly  crude)  example  is  to  partition  a  surveillance 
area  into  cells  and  assign  a  data  integration  center  to  each  cell.  In 
a  manner  similar  to  the  way  aviation  traffic  is  handled  by  air  traffic 
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control  centers*  sensor  data  for  one  geographic  cell  is  handled 
by  one  processing  center.  As  targets  leave  a  cell*  they  are  nhanded 
off"  to  another  center.  This  “cellular”  approach  has  the  advantage 
that  surveillance  data  flow  and  data  processing  is  more  evenly  distri¬ 
buted.  This  approach  also  has  advantages  in  flexibility.  Threat 
prioritization  and  resolution  can  be  handled  automatically  or  'hands-on* 
as  suits  the  command  for  a  region  and  the  resources  available  to  it. 

In  addition*  if  it  were  possible  that  all  centers  could*  at  their  option, 
tap  into  data  sources  for  other  cells  (by  fan-out  within  the  sensor 
data  network) ,  then  multiple  centers  could  focus  on  one  cell  for  higher 
system  performance  in  that  cell. 

The  accessibility  of  a  wide  range  of  sensor  data  to  a  large 
number  of  integration  centers  has  certain  implications  about  the 
sensor  data  network  associated  with  an  ocern  targeting  system.  An 
information  hierarchy  is  certainly  implied  in  that  there  is  a  wide 
mismatch  of  data  rates*  accuracy,  and  coverage  amongst  sensors. 
Existing  operational  sensor  systems  are  handled  for  the  most  part 
by  individual  baronies:  an  IR  satellite*  for  example,  downlinks  to  its 
own  earth  station  and  its  own  specialized  signal  processing  hardware. 
Therefore,  a  hierarchy  is  necessary  to  accommodate  existing  resources. 
One  way  to  provide  additional  crosslinking  between  sensors  and  com¬ 
putational  resources  could  be  via  communication  satellites  with  multiple 
spot  beam  antennas  for  both  ground  (or  air)  based  sensors  and 
users.  Intermediate  processing  for,  say,  an  IR  sensor  could  be 
accomplished  on  the  sensor  platform  or  at  a  ground  based  station. 
Whatever  results  of  the  processing  that  the  command  of  the  individual 
barony  wished  to  disseminate  could  be  crosslinked  to  other  sensor 
net  nodes. 

As  will  be  discussed  in  Section  8,  there  is  a  need  for 
data  bases  and  processing  at  the  force  element  level  as  well  as  at 
the  fleet  level.  Many  sources  of  information  may  be  fused  to 
form  a  fleet  level  picture.  Some  information*  such  as  from 
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wide  coverage  sensors  (satellite  based  staring  mosaics,  OTH  radar, 
etc.)  can  also  be  used  within  a  particular  tactical  region  (or  'cell') 
for  sensor  system  calibration.  An  aircraft  dropping  sonobuoys 
can  maintain  sonobuoy  position  relative  to  the  aircraft  by  local 
processing,  but  needs  at  least  inertial  navigation  input  to  maintain 
position  relative  to  other  units  and,  in  particular,  weapons  systems. 

e 

The  aircraft  knowledge  of  position  could  be  improved  if  its  local 
processor  could  tap  a  satellite  downlink  for  its  own  position.  A 
ground  based  acoustic  sensor  system  could  calibrate  its  own  coordinate 
system  with  the  fleet  coordinate  system  if  it  could  tap  the  location 
of  some  target  viewed  by  some  other  sensor  system  (OTH  radar,  etc.) 
as  well  as  its  own.  Exchange  of  data  and  data  processing  results  by 
cross-linked  satellites  (involving,  for  example,  FLTSATCOM)  and  by  a 
network  such  as  AUTODIN  II  could  open  up  a  host  of  distributed 
data  base  and  distributed  processing  possibilities  for  an  ocean  tactical 
targeting  system. 

There  are  several  difficult  areas  that  must  be  addressed  to 
solve  the  processing  problem: 

•  How  are  the  processing  problems  or  functions  to  be  parti¬ 
tioned  into  tasks  so  as  to  take  advantage  of  multiple 
processing  nodes? 

•  How  can  task  elements  be  automatically  identified  and 
assigned,  and  how  are  task  results  coordinated? 

•  If  the  allocation  of  tasks  is  to  be  distributed,  how  do 
various  processing  resources  negotiate  for  task  assignment? 

•  What  intermediate  representation  of  sensor  data  allows 
effective  utilization  by  a  wider  range  of  processing  functions 
within  the  targeting  system? 

•  In  what  ways  can  existing  sensor  baronies  be  interfaced  to 
a  network  in  such  a  way  that  their  output  is  more  generally 
useable  and  at  the  same  time  will  not  interfere  with  the  primary 
mission  of  the  sensor? 
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•  How  can  computer  hardware  and  software  systems  be 
designed  to  support  a  primary  function  specialization  as 
well  as  a  wide  range  of  temporary  specializations? 

•  If  it  is  the  case  that  it  is  not  practical  for  all  processing 
systems  to  have  access  to  all  data  (e.g.,  all  sensors  track 
file  from  section  2.0),  which  is  the  probable  case,  which 
sensor  systems  should  have  crosslinks  in  order  to  minimize 

duplicated  processing? 

I 

•  What  new  computer  technology  could  be  developed  to  handle 
the  extreme  response  time  requirements  for  threat  assessment/ 
threat  prioritization /weapon  allocation  problem  in  the  tactical 

^  theatre?  An  example  is  STARAN  type  associative  array 

processors  with  appropriately  configured  data  bases. 
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5.0  DATA  BASE  PROBLEMS 


When  considering  the  data  base  management  systems  to  be 
used  in  an  OT^  environment  with  many  sensors ,  it  is  essential  to 
approach  the  problem  at  the  same  conceptual  level  as  one  approaches 
the  rest  of  the  OT^  problem.  An  overall  system  design  should  be 
developed  that  will  integrate  the  various  data  types  that  will  be  used 
as  input.  The  design  should  take  into  account  the  relative  priorities 
of  different  data  bases  and  the  similarities  and  differences  in  their 
content  ^and  structure.  The  content  and  format  of  data  bases  are 
subject  to  change,  and  if  the  management  system  were  too  closely 
tied  to  specific  data  base  configurations,  any  change  in  configuration 
would  require  extensive  revisions  to  the  data  base  management  system 
as  well. 

In  the  lower  levels  of  the  OT^  hierarchy,  there  are  several 
requirements  for  the  handling  of  data  in  its  raw  and  near-raw  state. 
TARF  data  must  be  filed  in  its  raw  form  upon  receipt,  but  error 
correcting  and  field  validation  processing  might  also  imply  a  parallel 
TARF  data  base  with  'corrected'  entries.  Time  series  datal  weather 
satellite  data,  and  OTH  radar  data  all  require  data  bases  designed 
for  the  particular  data  rate  characteristics,  burstiness,  and  access 
methods  needed  for  the  data  involved.  Because  the  data  bandwidth 
of  sensor  data  is  larger  near  the  sensor,  more  attention  must  be  paid 
to  sensor  data  base  performance.  T 

Since  a  data  base  maps  data  flow  from  a  source  process  to  a 
sink  process  (as  well  as  acting  as  a  repository),  the  access  pattern 
or  'key1  used  to  store  the  data  is  often  different  from  the  retrieval 
pattern.  For  this  reason,  access  to  a  given  data  base  should  not  be 
by  a  single  key  or  by  a  hierarchically-ordered  set  of  keys.  Such  an 
access  plan  assumes  that  requests  to  the  file  for  information  will 
always  be  of  the  same  kind  (i.e.,  by  STAR  target  number),  and  is 
incapable  of  efficient  retrieval  by  lower-ordered  keys.  A  number  of 
data  access  methods  exist  which  handle  multiple  keys  efficiently. 


These  methods  can  be  tuned  to  give  optimal  retrieval  time  for  the 
most  commonly  used  keys  without  sacrificing  general  efficiency  no 
matter  what  key  or  combination  of  keys  is  specified. 
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One  such  access  method  is  partial  match  retrieval  using  hash¬ 
ing  and  trie  algorithms.  This  is  an  associative  retrieval  method  with 
proven  upper  bounds  on  the  number  of  disk  accesses  required  to 
retrieve  information  specified  with  varying  numbers  of  keys. 

In  all  hashing  algorithms,  a  function  is  applied  to  the  original 
set  of  key  values  to  create  a  new  value  which  is  the  address  in  the 
file  where  the  item  will  be  stored.  To  retrieve  the  item,  the  key  is 
hashed  with  the  same  function  to  produce  the  item's  address.  The 
unique  aspect  of  partial  match  algorithms  is  that  they  can  retrieve 
data  items  which  have  specific  values  for  some  keys  and  unspecified 
or  wild  card  values  for  other  keys. 

A  system  of  this  kind  has  very  desirable  properties  for  a  par- 

2 

ticular  application  like  an  OT  test  bed.  Insertion  of  new  data  is  very 
fast.  The  system  permits  access  to  data  bases  on  an  assortment  of 
keys  and  facilitates  retrieval  of  groups  of  data  items.  For  example, 
in  dealing  with  input  data,  an  operator  might  wish  to  confine  atten¬ 
tion  to  peaks  above  a  stated  coherence  level  whose  time  of  receipt 
was  within  a  certain  time  slot,  or  peaks  from  a  given  pair  of  sensors, 
or  peaks  from  a  given  geographical  part  of  the  ocean  that  satisfy 
various  fequency  and  water  time  constraints.  If  this  data  is  organized 
by  a  single  key  or  by  a  hierarchically-ordered  set  of  keys,  few  queries 
of  this  kind  can  be  handled  efficiently,  it  at  all.  A  multiple- key , 
partial-match  hashing  algorithm  handles  queries  of  this  nature  easily. 
When  the  software  receives  a  retrieval  request,  specified  key  values 
are  left  intact  and  values  for  unspecified  keys  are  treated  as  wild 
cards . 

Another  multiple  key  access  method  employs  a  Normal  Multi¬ 
plication  Table  as  an  index  to  the  data  base.  This  method  is  well- 
suited  to  data  bases  where  the  number  of  data  items  is  large  relative 
to  the  number  of  distinct  key  values. 
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As  was  the  case  with  trie  algorithms,  system  performance  can 

be  tuned  by  ranking  the  keys  according  to  expected  fequency  of 

use.  singly  or  in  combination.  Perhaps  the  greatest  difference 

between  the  two  methods  is  that  the  Normal  Multiplication  Table  is 

constructed  from  the  contents  of  a  specific  data  base  while  the  trie 

algorithm  is  designed  to  handle  any  data  set  that  contains  items  that 

can  be  described  by  the  current  set  of  keys.  This  difference  has 

2 

important  implications  for  an  OT  test  bed.  The  Normal  Multiplication 
Table  method  would  work  well  for  artificial  input  data  or  for  input  data 
that  arrived  in  discrete  blocks,  as  is  the  case  with  the  "epochs"  cur¬ 
rently  being  used  in  experiments  at  the  ARC.  The  algorithms  that 
build  the  table  are  relatively  simple  and  fast,  and  once  the  table  is 
built,  it  would  provide  extremely  fast  information  retrieval. 

Because  they  are  attribute-based,  both  of  these  methods  adapt 
easily  to  changes  in  file  format  or  content.  They  require  far  less 
core  space  than  the  standard  inverted  indices  or  multilist  directories. 

If  proper  data  definition  is  done  at  the  outset,  these  methods  can  be 
used  across  files  and  therefore  can  link  together  data  from  different 
input  data  bases. 

Information  dealt  with  at  higher  levels  in  an  OT^  system  has 

been  extracted  from  data  at  the  sensor  by  various  processing  tech- 
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niques.  Higher  level  OT  data  bases  will  contain  such  information  as 
target  identification,  position,  course,  and  speed.  Target  identification 
may  include  hull  type,  known  weaponry,  intent,  and  severity  of 
threat  if  an  enemy.  In  addition  to  representing  the  current  state  of 
the  surveillance  area,  a  high  level  data  base  may  contain  a  target 
history  of  past  movements  and  events,  and  maintain  a  list  of  sensors 
contributing  to  a  piece  of  high  level  information. 

A  commander  may  need  to  relate  the  positions  of  certain  fleet 
elements  to  the  sea  condition  for  that  particular  area.  However,  the 
location  of  a  data  base  containing  weather  may  be  physically  separated 
from  the  command  center  data  base.  Therefore,  issues  concerning 
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distributed  data  bases  need  to  be  addressed.  Various  fusion  centers 
will  maintain  data  bases  relevant  to  their  own  sensor  suite.  However, 
decisions  made  at  a  central  location  require  that  all  these  data  bases 
be  accessible  in  a  uniform  manner.  In  addition,  procedures  will  be 
required  for  the  event  of  the  failure  of  a  remote  data  base. 

The  wide  range  of  types  and  formats  of  data  required  at  the 
fleet  level  may  argue  for  a  relational  data  base  design.  Ultimately, 
one  would  wish  to  automate  fairly  complex  queries  of  the  data  base 
to  identify  particularly  important  situations  without  human  intervention. 
The  structure  of  a  relational  data  base  is  well  suited  to  the  associative 
searches  that  are  required  to  handle  queries  of  this  kind.  Because 
the  information  is  organized  by  its  attributes  and  does  not  have  a 
structure  imposed  on  it,  it  can  be  efficiently  searched  for  records 
which  match  a  combination  of  key  values  that  characterize  a  particular 
threat  or  other  tactical  situation.  In  instances  where  speed  is  critical, 
special  hardware,  such  as  fast  array  processors,  can  be  used. 

2 

For  the  highest  levels  in  an  OT  system,  relational  /distributed 
data  base  technologies  exist  that  can  handle  the  sorts  of  information 
required  by  a  commander  at  the  fleet  level.  Some  appropriate  work 
with  human  interface  technology  and  relational  data  bases  has  been 
accomplished  through  the  ACCAT  project.  However,  serious  issues 
of  system  performance  have  to  be  addressed.  Also,  it  is  likely  that 
several  data  base  access  techniques,  some  of  which, weje  discussed 
above,  would  be  more  suited  for  data  at  the  local  tactical  level 
because  of  response  time  constraints.  The  nature  of  the  data  being 
dealt  with  at  levels  closer  to  the  sensor  is  also  quite  different.  The 
point  is  that  though  the  interfaces  to  different  data  base  systems 
should  be  similar,  the  interfaces  should  be  independent  of  the  access 
methods  used.  The  access  methods  should  depend  on  the  type  of 
data  handled,  types  of  access,  and  response  time  requirements.  A 
uniform  interface,  on  the  other  hand,  will  lend  itself  to  uniform 
representation  of  sensor  information  and  thus  more  flexible  utilization 
of  data  throughout  the  targeting  system. 
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Thus*  several  problem  areas  that  must  be  addressed  to  properly 

2 

use  data  base  technology  in  an  OT  system  are: 

•  Identification  of  what  types  of  data  bases  are  actually 
required  for  the  OT^  problem.  Examples  range  from 
acoustic  data  bases  at  an  acoustic  data  processing  center 
to  data  bases  containing  hull  number  and  ship  captain 
for  a  particular  ship  at  the  fleet  level. 

•  What  data  base  access  techniques  are  optimal  for  data  at 
various  levels  in  the  hierarchy? 

•  Can  specialized  hardware  be  applied  to  certain  time-critical 
applications?  Cellular-logic  devices*  head-per- track  rela¬ 
tional  data  base  machines*  associative  array  processors*  and 
mass  storage  devices  utilizing  content  addressable  blocks 
are  examples. 

•  At  what  level  in  the  system  are  distributed  data  bases 
appropriate?  What  types  of  information  should  he  distri¬ 
buted  and  what  types  tend  to  be  centralized?  How  can 
consistency  between  such  data  bases  be  maintained?  How 

is  the  problem  of  failure  of  remote  sections  of  the  data  base 
handled  (fail  soft)? 

•  For  time- dependent  data*  what  techniques  may  be  used  to 
identify  and  delete  unimportant  data  from  the  data  base 
in  order  to  avoid  data  base  overflow? 
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6.0  MEASUREMENT  FUSION 


In  considering  the  fusion  problem  embedded  in  the  OT^  problem, 
one  must  give  paramount  importance  to  the  data  available  and  its  utility. 
In  concept  formulation  one  should  not  be  tied  too  closely  to  any  given 
set  of  sensors  or  sensor  system  but  rather  to  concepts  and  potential 
suites  of  sensors.  There  are  a  large  number  of  possible  sensors 
whose  data  could  be  potentially  useful  in  the  fusion  center.  Types 
of  sensors  that  could  be  considered  are: 

•  Acoustic ; 

•  Electro-magnetic; 

•  Optical ; 

•  Infra-red; 

•  Thermal;  and 

•  Human  sources. 

In  most  cases  both  passive  and  active  sensors  could  be  con¬ 
sidered  as  well  as  both  fixed  and  mobile  locations  and  many  possible 
sensor  placements  (satellite,  land,  air  or  sea). 

With  this  plethora  of  possible  data  and  data  types  many  ques¬ 
tions  and  problems  arise.  In  the  simple  case  of  detecting  one  signal 

in  noise  or  of  tracking  a  given  target  (with  data  for  that  target) 

*  ^ 

the  addition  of  more  data  either  helps  or  does  not  hurt  the  processor. 
This  is  not  true  of  the  data  association  problem,  particularly  with 
limited  communication,  computational  or  human  resources.  It  is  possible 
in  the  multi-target/sensor  arena  to  turn  a  situation,  where  one  could 
detect,  track  and  classify  a  number  of  targets,  into  one  where  this 
is  impossible.  This  could  be  made  to  occur  for  any  given  system  by 
the  addition  of  more  clutter,  more  real  targets,  or  even  more  good 
data  from  existing  targets.  This  could  occur  because  of  the  clogging 
of  limited  communication  channels,  the  creation  of  too  many  possible 
or  feasible  combinations  of  data  in  a  multi-target  tracker  or  saturation 
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of  computational  resources.  In  this  somewhat  novel  situation  it 
becomes  very  important  to  define  the  quality  measure  that  one  will 
use  to  evaluate  the  eventual  system  output.  It  must  be  noted  that 
there  could  be  multiple  (possibly  conflicting)  requirements  on  this 
system  requiring  a  multi- dimensional  quality  measure. 

Since  one  must  finally  operate  the  data  association,  detection, 
classification,  localization  and  tracking  systems  on  a  set  of  data  from 
a  set  of  sensor  types,  one  can  think  of  each  possible  selection  of  data 
as  one  set  of  possible  or  available  data.  In  evaluating  a  given  set  of 
some  allowable  communication /computer  resource  suites,  one  could 
choose  a  number  of  quality  measures.  One  could  perhaps,  from  a 
technical  point  of  view,  define  the  information  content  of  a  set  of 
measurements — some  kind  of  entropy  measure  on  an  absolute  basis. 

A  measure  of  a  set's  ability  to  give  high  resolution  target  tracks 
(most  precise  tracks)  could  be  one  such  quality  measure.  Another 
would  be  a  measure  of  the  ability  of  a  set  to  allow  separation  of 
different  targets  or  different  classes  of  targets.  In  this  case  one 
could  look  at  the  features  that  one  wanted  to  separate,  define  a 
distance  measure  and  use  pattern  recognition  techniques.  One  might 
opt  for  the  set  that  would  allow  the  largest  number  of  targets 
(possibly  of  a  single  class)  to  be  detected  while  minimizing  the 
probability  of  false  alarm. 

It  should  be  clear  that  if  a  system  used  its  resources  to  track, 

*  *  » 

detect  or  classify  loud  targets  it  might  not  have  the  resources  required 
to  track  quieter  or  other  particularly  important  targets.  The  measure¬ 
ment  sets  allowing  best  detection  might  not  even  be  observable  sets 
for  track  reconstruction  or  target  classification  or  separation  (data 
association) . 

Other  means  of  measurement  quality  could  include  time  late¬ 
ness,  independence  of  other  measurements  in  the  set,  and  correlatability 
with  the  measurements  from  different  systems.  Even  if  some  measure 
of  goodness  data  from  one  sensor  system  is  judged  to  be  the  best. 
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additional  data  from  that  system  on  a  given  target  would  not  be 
independent.  Thus,  it  might  be  better  to  include  data  (even  of  a 
lower  quality)  for  a  different  sensor  system  to  confirm  the  existence 
of  a  given  target. 

Major  research  problems  exist  in  the  data  fusion  area.  Some 
of  these  are  data  independent  but  many  others  depend  on  the  details 
of  the  sensor  system,  the  statistical  characteristic  of  the  measure¬ 
ment  systems  and  even  cost  of  data  transmission,  computation  and 
other  factors,  allowing  a  reasonable  comparison  of  resources  on  a  cost 
or  uniform  basis. 
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7.0  MULTI-SENSOR /MULTI-TARGET  TRACKING,  DETECTION 

AND  CLASSIFICATION  PROBLEMS 

The  multi-sensor  /multi-target  tracking /detection /classification 
problem  area  is  a  formidable  one.  For  a  given  surveillance  area,  there 
may  be  a  large  number  of  surface /subsurface  ships,  aircraft,  and  other 
acoustical  signal  sources  in  close  proximity  to  one  another  such  that 
one  is  confronted  with  a  dense  multi-target  environment. 

To  illustrate  the  complexity  of  the  problem,  one  could  first 
consider  the  case  where  there  is  one  sensor  system  that  is  receiving 
signals  from  several  targets  as  illustrated  in  Figure  6  below,  where 
each  of  the  targets  could  potentially  represent  aircraft,  surface/ 
subsurface  ships,  etc.  The  sensed  signal  might  even  be  caused  by 
sea- induced  noise. 


SENSOR 

SYSTEM 


Figure  6.  The  multi- target  problem. 
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This  problem  is  defined  as  the  single-sensor,  multi-target 
problem.  When  the  set  of  targets  (A,B,C,D)  are  operating  in  close 
quarters  and/or  when  the  sensed  signals  have  a  low  signal-to-noise 
ratio,  a  fundamental  problem  is  that  of  data  association  (return-to- 
track).  That  is,  for  each  incoming  measurement,  how  do  we  attach 
it  to  the  correct  target  that  produced  this  measurement,  or  is  the 
measurement  caused  solely  by  noise  (false  alarm  or  clutter)  as  in 
box  E  in  the  above  diagram.  This  is  fundamentally  a  target 
detection  problem.  There  have  been  a  number  of  approaches  used 
in  the  past  in  solving  this  problem,  and  they  have  met  with  varying 
degrees  of  success.  They  range  in  scope  from  very  simple  procedures 
such  as  manual  selection  of  data  by  a  human  operator  to  more  compu¬ 
tationally  sophisticated  methods  that  can  automatically  handle  the  data 
association  problem  without  operator  intervention. 

The  manual  procedures  involve  careful  preprocessing  of  the 
data  by  the  operator,  and  although  laborious  and  time-consuming, 
have  met  with  some  success  on  real  data.  They  work  well  when 
the  targets  are  widely  separated  and  the  signals  have  a  high  signal- 
to-noise  ratio,  but  in  a  "dense"  environment  and/or  low  signal-to- 
noise  ratios,  these  approaches  are  not  very  suitable.  It  is  the  latter 
set  of  automated  algorithms  that  fare  better  under  this  type  of  scenario, 
and  lately  have  been  shown  to  exhibit  quite  good  tracking  performance 
on  both  real  and  simulated  data.  Graphic  examples  of  this  were  found 
during  the  SASE  experiment. 

For  the  most  part,  all  of  the  automated  methods  employ  extended 
Kalman  filtering  of  one  form  or  another  to  perform  data  association. 

This  is  done  using  statistics  generated  by  the  filter  and  subsequently 
comparing  them  to  a  set  of  specified  thresholds.  In  addition  to  the  data 
association  task,  the  filter  is  structured  to  do  the  target  tracking  in 
parallel  so  as  to  provide  target  estimates  in  a  preferred  coordinate 
system  (for  instance,  position,  speed  and  course). 


In  addition  to  the  problems  already  noted  in  the  above  dis¬ 
cussion,  one  must  now  face  the  issue  of  correlating  the  measurements 
from  different  sensor  systems  for  the  same  given  target.  This  is  a 
multi-sensor  measurement  fusion  problem  and  a  few  solution  approaches 
to  it  have  already  been  suggested;  namely,  independent  multi-target 
tracking  for  each  sensor  system,  and  then  the  combination  of  the 
individual  sensor  tracks  into  a  resultant  set  of  target  tracks. 

To  illustrate  the  complexities  involved  in  going  to  a  more 
elegant  approach  that  would  be  an  improvement  to  the  above,  it 
might  be  helpful  to  pose  a  specific  approach  as  an  example.  This 
approach  would  have  the  tendency  to  increase  the  dimensionality 
of  the  observation  space  in  a  sequential  fashion  thereby  increasing 
the  information  content  of  the  data  set.  The  key  features  can  best 
be  described  with  the  aid  of  the  block  diagram  in  Figure  8.  Essen¬ 
tially  it  involves  implementing  a  multi-target  tracking  filter  to  per¬ 
form  the  detection,  classification,  and  tracking  functions  for  each 
of  the  individual  sensor  systems  A  to  2,  respectively.  The  tracks 
generated  by  each  of  the  filters  are  then  combined  into  a  composite 
set  of  tracks.  If  a  given  target  is  being  sensed  by  more  than  one 
sensor  system,  the  filter  for  one  sensor  system  should  generate  a 
target  track  that  is  "reasonably  close"  to  a  target  track  generated 
by  the  filter  for  another  sensor  system.  The  major  objective  here 
is  to  correlate  those  target  tracks  of  different  sensor  systems  that 
correspond  to  the  same  target.  This  is  the  function  of  the  track 
clustering  block  shown  in  Figure  8. 

A  variety  of  clustering  approaches  can  be  used  to  identify 
tracks  belonging  to  a  common  target.  They  range  in  complexity 
from  simple  procedures  such  as  a  least  squares  fit  to  more  complex 
procedures  that  use  the  confidences  associated  with  a  set  of  target 
tracks  to  "weight"  the  results  of  the  clustering  approach.  This 
latter  approach  can  be  implemented  recursively  using  an  extended 
Kalman  filter  which  has  as  input  measurements  the  target  track 
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When  using  the  filter  for  data  association  and  target  tracking, 
care  must  be  taken  in  selecting  a  filter  model  that  accurately  reflects 
the  physics  of  the  actual  target  motion.  The  usual  assumption  is 
constant  rhumb  line  or  great  circle  motion.  This  is  normally  adequate 
unless  the  target  starts  maneuvering.  When  this  occurs,  the  filter 
could  misinterpret  measurements  from  a  maneuvering  target  for  clutter 
when  performing  the  data  association  task,  and  not  use  the  measurement 
at  all  in  forming  a  target  track.  For  abrupt  maneuvers,  solutions  to 
this  problem  can  normally  be  achieved  by  artificially  increasing  the 
filter's  process  noise  covariance  matrix.  However,  for  more  subtle 
maneuvers  an  obviously  more  sophisticated  algorithm  must  be  applied. 
Several  innovative  approaches  have  been  proposed  recently  to  differentiate 
between  maneuvers  of  this  sort  and  clutter.  One  of  the  more  promis¬ 
ing  of  these  appears  to  be  the  accumulation  of  frequency  sensitive 
measurement  residuals  as  a  test  statistic  for  maneuver  detection. 

The  above  discussion  has  thus  far  only  considered  a  single 
sensor  system.  A  much  more  complex  problem  is  introduced  by  con¬ 
sidering  the  multi-sensor  system,  multi-target  problem  depicted  in 


Figure  7.  The  multi-sensor /multi-target  problem. 


estimates  and  covariance  matrices  generated  by  the  filters  for  each 
of  the  sensor  systems  A  to  Z,  respectively.  This  comprises  the 
inner  feedback  loop  of  the  track  clustering  algorithm  shown  in 
Figure  8. 

Having  obtained  a  set  of  clustered  tracks,  they  can  then  be 
transformed  back  into  a  set  of  estimates  in  measurement  coordinate 
units  for  each  of  the  sensor  systems.  Association  of  the  actual 
measurements  with  the  measurement  estimates  above  is  obtained  by 
identifying  those  measurements  lying  "close11  to  the  estimates  for  a 
given  target  track  from  the  cluster  for  each  of  the  sensor  systems. 

The  measurements  for  the  given  target  track  are  then  cor* 
related,  i.e. ,  they  are  identified  as  having  come  from  the  same  target, 
and  are  combined  into  an  increased  dimension  observation  space  having 
more  information  content  than  the  measurements  for  each  of  the  sensor 
systems  by  themselves.  A  filter  can  then  be  designed  to  process  this 
multi-sensor  data  to  produce  a  more  accurate  target  track.  This 
procedure  continues  until  all  member  tracks  of  the  target  track  cluster 
have  been  processed. 

At  this  point,  the  outer  feedback  loop  on  the  track  clustering 
algorithm  can  be  closed  to  produce  a  refined  set  of  target  tracks. 

This  is  done  by  repeating  the  track  clustering,  measurement  data 
association,  and  multi-sensor  filtering.  The  only  difference  here 
is  that  the  track  clustering  is  initiated  on  a  set  of  composite  tracks 
that  were  generated  by  the  multi- sensor  filter.  It  is  easy  to  see 
how  this  can  be  put  in  an  iterative  mode,  i.e.,  continue  to  feed  back 
the  resultant  target  tracks  into  the  track  clustering  algorithm,  until 
no  further  change  in  the  resultant  target  tracks  is  noted. 

Again,  this  approach  has  been  intended  as  an  example  and 
not  as  a  specific  solution  to  the  multi- sensor  measurement  fusion 
problem.  The  purpose  here  is  to  illustrate  the  intricacies  associated 
with  any  approach  that  would  extract  all  of  the  information  embedded 
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in  a  set  of  measurements  from  different  sensor  systems,  and  use  it 
in  an  optimal  fashion  to  enhance  our  knowledge  in  the  areas  of 
detection,  classification,  localization  and  tracking  of  target  threats 
within  the  surveillance  area. 


I 


I 


I 


Cr 

l 

I 


» 


66 


ORINCON 


I 


4 


> 


I 


I 


» 


» 


» 


I 


s 


» 


» 


8.0  FORCE  AND  WEAPON  ALLOCATION 

The  ultimate  goal  of  an  Ocean  Tactical  Targeting  System  is 
to  bring  the  proper  weapons  to  bear  on  the  targets  in  order  to 
neutralize  them  while  minimizing  the  losses  among  our  own  forces. 
There  are  three  issues  which  need  to  be  addressed  in  this  area: 

e  Threat  prioritization ; 

•  Friendly  force/target  assignments;  and 

•  Real-time  battle  management. 

As  part  of  the  threat  prioritization  function,  the  data  base 
management  system  will  provide  the  information  to  allow  priorities  to 
be  assigned  to  each  of  the  threat.  This  prioritization  will  be 
assigned  based  on: 

•  Threat  type ; 

•  Location;  and 

•  Intent  (e.g.,  direction  of  movement). 

For  example,  a  high  priority  will  be  given  to  an  enemy  ship  carrying 
cruise  missiles  near  a  friendly  land  mass,  whereas  a  low  priority 
might  be  assigned  to  an  enemy  ship  carrying  only  surface-to-air 
missiles  located  far  from  any  friendly  ships  or  land  masses. 

In  the  actual  targeting /weapon  allocation  problem  there  are 
hierarchical  levels  of  command  and  control  which  must  be  considered 
for  effective  use  of  the  available  weapon  resources.  The  top  level 
of  this  structure  (see  Figure  9)  will  allocate  groups  of  threats  to 
friendly  forces  based  on  the  threat  priorities.  As  time  evolves,  the 
threats  will  move,  their  assigned  priorities  will  change,  and  the 
force /target  allocations  will  be  modified.  These  changes  will  take 
place  relatively  slowly  so  that  these  assignments  can  be  made  in  an 
optimal  manner. 
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The  force  level  battle  management  function  takes  place  at  a 
lower  level.  This  includes  real-time  ship  or  weapon  assignment  to 
threats,  and  real-time  reallocation  of  forces /weapons  based  on  battle 
assessment  data.  This  level  is  further  complicated  by  the  need  to 
counter  enemy  attacks,  as  well  as  to  destroy  the  assigned  threats. 
Because  of  the  real-time  requirement  of  this  battle  management 
function,  some  degree  of  optimality  in  weapon  allocation  must  be 
sacrificed  to  achieve  the  real  time  capability. 

Such  a  two-level  command  and  control  structure  allows 
optimality  of  force /threat  allocations  to  be  emphasized  because  of  the 
relatively  slow  dynamics  at  this  upper  level  while  delegating  to  a 
lower  level  the  real-time  battle  management  function. 

This  is  a  complicated  resource  allocation  problem  and  there 
are  a. number  of  specific  issues  which  need  to  be  addressed: 

•  What  is  the  best  division  of  functions  between  the 
two  levels? 

•  How  detailed  must  the  weapon  allocation  be,  e.g.,  ship  level 
or  weapon  level? 

•  How  are  the  offensive  and  defensive  issues  of  battle  man¬ 
agement  optimally  addressed? 

•  How  can  these  approaches  be  automated,  and  can  the 
concepts  of  artificial  intelligence  be  used? 

This  might  allow  real-time  planning  to  improve  battle  management 
at  a  lower  level,  freeing  the  major  human  decision-makers  to  make 
major  policy  decisions  or  system  overrides. 
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9.0  FINAL  TARGETING 


When  the  long  range  surveillance  system  has  detected,  tracked 
and  classified  a  target  to  the  best  of  its  ability,  one  generally  will 
still  not  know  the  target  position  well  enough  to  launch  a  weapon 
system  against  it.  It  may  also  be  politically  unacceptable  to  launch 
an  over-the-horizon  weapon  without  a  final  classification  system 
that  will  confirm  the  target  classification  made  at  long  range.  There 
are  many  unsolved  research  problems  inherent  in  end  point  targeting 
and  classification.  The  first  targeting  could  be  done  using  a  deploy¬ 
able  sensor  (perhaps  a  sonobuoy  field)  in  conjunction  with  one  or 
more  aircraft  and  a  weapon  system. 

A  strike  system  designed  to  destroy  a  submarine  must  obtain 
positional  information  about  the  submarine  relative  to  the  strike  air¬ 
craft.  In  order  to  ascertain  this  information,  a  sonobuoy  field  may 
be  used.  The  field  must  contain  enough  buoys  to  obtain  sufficient 
observables  for  determination  of  the  necessary  position  and  velocity 
information.  The  optimum  number  is  scenario  dependent.  The  desig¬ 
nation  aircraft  which  accomplishes  the  target  function  must  derive 
the  information  from  the  observables  and  transfer  this  information 
to  the  strike  aircraft  which  must,  in  turn,  establish  an  attack  posi¬ 
tion  and  launch  its  weapon.  The  type  of  weapon  guidance  used  is 
also  scenario  dependent.  The  weapon  must  enter  the  water  near 
enough  to  the  submarine  to  allow  the  terminal  sensor  to  find  and 
lock  onto  it.  A  stealthy  attack  is  required  in  order  to  negate  the 
target's  ability  to  employ  countermeasures  or  counterthreats  and 
to  deny  performance  of  its  main  function.  Thus,  a  detailed  analysis 
of  the  scenarios,  function  and  operational  capability  of  such  a  sonobuoy 
field  target  strike  system  must  be  performed.  From  earlier  ORINCON 
studies  performed  for  NADC,  it  appears  that  there  may  be  fundamental 
observability  problems  which  must  be  investigated.  Problem  areas 
include  the  strike  scenario,  the  prediction  of  target  position  and 
velocity,  the  weapon  guidance  modes,  transfer  alignment  and  the 
effect  of  the  terminal  sensor. 
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The  problem  of  utilization  of  a  sonobuoy  field  to  obtain  target¬ 
ing  information  and  to  perform  a  strike  mission  against  a  submarine 
using  guided  weapons  requires  a  detailed  consideration  of  system 
options  and  their  accuracies.  Targeting  utilizes  the  sonobuoy  refer¬ 
ence  system  to  obtain  the  position  and  velocity  of  the  submarine 
relative  to  the  aircraft.  This  includes  the  aircraft  Doppler  and 
inertial  navigation  systems,  and  an  azimuthal  measurement.  The 
interferometric  angular  determination  requires  an  antenna  pair  base¬ 
line  on  the  aircraft  and  a  transmitter  on  the  sonobuoy.  The  targeting 
of  the  submarine  has  uncertainties  involved  due  to  the  aircraft  inertial 
system  and  interferometer  errors.  These  cause  a  degradation  of  the 
submarine  location  and  velocity  relative  to  the  designation  aircraft. 

It  has  been  found  that  with  the  use  of  sonobuoys  this  designation 
problem  is  highly  sensitive  to  aircraft  maneuvers. 

When  the  designation  aircraft  is  not  the  strike  aircraft,  further 
problems  occur.  In  particular,  the  strike  aircraft  position  relative  to 
inertial  space  is  in  error  due  to  its  inertial  navigation  system  errors. 
Thus,  since  the  designation  aircraft  position  relative  to  inertial  space 
also  contains  uncertainties,  the  location  of  the  strike  aicraft  relative 
to  the  designation  aircraft  is  in  error.  This  implies  an  additional 
uncertainty  in  the  location  and  velocity  of  the  submarine  relative  to 
the  strike  aircraft. 

The  performance  of  the  strike  system  may  be  affected  in  various 
ways.  Of  course,  the  basic  intent  and,  thus,  the  ultimate  performance 
measure  is  that  of  destroying  the  submarine.  In  order  to  accomplish 
this,  the  weapon  must  be  placed  within  a  certain  volume  about  the 
submarine  in  order  to  assure  that  the  terminal  sensor  can  activate 
and  lock  onto  the  target.  Failure  to  place  the  weapon  within  this 
boundary  will  result  in  an  unsuccessful  strike,  and  is  due  to  the 
system  errors  discussed.  These  errors  are  dependent  upon  strike 
parameters  such  as  launch  range,  type  of  guidance  laws  and  systems, 
number  of  sonobuoys  in  the  field,  proximity  of  strike  aircraft  to 
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designation  aircraft  and  ability  of  the  designation  aircraft  to  maneuver 
to  obtain  the  proper  observability. 

Another  performance  measure  is  that  of  the  ability  of  the 
aircraft  to  survive  any  counterthreat.  This  may  require  a  "launch 
and  run"  scenario  rather  than  an  "on-station  until  destroy"  scenario. 
The  ability  of  the  strike  system  to  fire  at  long  range  is  certainly 
important  in  the  counterthreat  environment.  Other  measures  include 
the  ability  to  distinguish  real  from  imaginary  threats ,  to  acquire, 
designate  and  strike  quickly  and  reliably. 
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